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ABSTRACT 


The tidal perturbations in the orbital Inclinations of the BE-C, GEOS-I, 
and GEOS-H satellites are analyzed. Effective tidal Love numbers 
and phase angles for the , M 2 , Kj, Pj , and Sj tides are recov- 

ered. The effective tidal phase angles tend to be on the order of a few 
degrees. The effective tidal Love numbers are generally less than the 
solid earth Love number k, of about 0.30. This supports the contention 
of Lambeck, et al. (1974) that the ocean tides give an apparent depres- 
sion of the solid earth Love number. Ocean tide amplitudes and phases 
are calculated for the above tides assuming kj * 0.30 and the solid 
earth lag angle » 0. The results show good agreement with 
Felsentreger, et al. (1976) on GEOS-I but not on GEOS— II, The 
effective Love number and phase angle are poorly determined, but give 
a lunar acceleration of r^29±15 arc sec/(100yr)^, an energy dissi- 
pation of -3.6 ±1.8 X 10^^ erg/sec, and a tidal function time scale 
of 1.4 X 10^ yr when averaged over all three satellites. This is in 
fair agreement with current estimates. 
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TIDAL PARAMETERS DERIVED FROM THE PERTURBATIONS IN 


THE ORBITAL INCLINATIONS OF THE BE-C, GEOS-I, 
AND GEOS-II SATELLITES 


INTRODUCTION 

The sun and moon raise tides on the earth; the displaced mass perturbs 
the orbits of earth satellites. These perturbations may be analyzed in order to 
learn about the earth's tides. 

Early analyses of the tidal perturbations performed by Kozai (1965, 1968), 
Newton (1965, 1968), Anderle (1971), Smith, et al. (1973), and Douglas, et al. 

(1974) were aimed at recovering the Love number k, and lag angle of the 
solid earth tide. The resulting Love numbers were found to be generally 
smaller than the well-known value for k^ of about G.29 to 0.31 obtained from 
seismic studies (e.g., Longman, 1966). The discrepancy has been ascribed by 
Lambeck (1974) and Lambeck et al. (197 4) to the corruption of the solid earth 
tidal signal by the ocean tides (and to a lesser extent the atmospheric tides), 
which Kaula (1962) showed might be significant. The ocean tides give the appear- 
ance of Love numbers and lag angles which vary from constituent to constituent. 
Accordingly, some satellite data have been re-analyzed (Lambeck, et al. 1974; 
Felsentreger, et al. 1976) in an effort to learn about the solid earth, ocean, and 
atmospheric tides. 
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We continue with such a re-analysis here by examining the tidal perturba- 
tions in the orbital incUnation of three different earth satellites: BE-C (Beacon 
Explorer C), GEOS-I, and GEOS-II. Table 1 gives pertinent information about 
the satelUtes. Inclination perturbations are analyzed because they give a rela- 
tively clean separation of tidal effects from other disturbing forces. 

THE DATA 

The BE-C data consists of 36 laser observations, spanning a 501 day period 
beginning in July 1970. The amplitude of the variation in inclination is about 1.0 
arc seconds (40 meters). For further details, see Smith, et al. (1973). Tlie 
GEOS-I consists of 142 TRANET Doppler observations spaced over 626 days, 
beginning in February 1966. The amplitude of the tidal signal is about 1.2 arc 
V, second (36 meters). The GEOS-H data is composed of 113 precision reduced 
camera and TRANET Doppler observations taken over a 651 day period begin- 
ning in March 1968. For further details, see Felsentreger, et al. (1976). 

In each case the well-known disturbing forces were modelled and removed, 
so that the remaining perturbations are presumably due almost entirely to the 
tides, . 

ANALYTICAL EXPRESSION FOR THE INCLINATION VARIATION 

We now derive an analytical expression for the total tidal perturbation in 
the inclination of the orbit of an earth satellite by modifying the equation for the 
solid earth tidal potential. 
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The solid earth tidal potential expressed in orbital elements is given by 


(Kaula, 1964) 


^ftmpq ^6mpq ^ ^^nihj 

ii.i 


• cos 


^2mpq 




( 1 ) 


where 


:| 

Km 


= Gm* 


(C - m) ! 
(K+ m)! 



Ftn.pO*)Ctp,<«*) 


and 


'^Empq ~ “ 2p)0J* + (G - 2p + q)M* + 


The starred (*) quantities refer to the disturbing body (sun or moon) and the 
unstarred refer to the satellite. The sign in front of the tidal phase angle 
had been minus iuKaula's original formulation; Lambeok, et al. (1974) changed 
the sign to plus, which makes the important phase angles positive quantities in 
the ease of a frictionally-delayed solid earth tidal bulge. We follow the conven- 
tion of Lambeck, et al. here. 

We wish to find the total tidal potential Tg^^^ caused by the solid earth, 
oceans, and atmosphere, since this is the potential felt by a satellite. We limit 
the discussion to 2 =* 2, since the second degree terms dominate the motion of 
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the satellite and the evolution of the earth-moon system. (For a discussion of 
fourth degree terms, see below and Appendix 2.) 

In finding the total tidal potential of degree 2 we make the following assump- 
tion: the amplitude of the potential and the phase angle become frequency- 
dependent. This means that the solid earth Love number k, in eq. (1) is replaced 
by the effective Love number and the solid earth phase angle is 

replaced by the effective phase angle . The adjective ’’effective" is neces- 

sary since we now refer to the tidal response of the entire earth, and not just 
the solid earth. This procedure takes care of any possible frequency dependence 
of the solid earth Love number, the response of the oceans, and the atmosphere 
tides. The letters "d" and” 6" are chosen to avoid possible confusion with the 
solid earth tides. 

The total tidal potential is now to be substituted in Lagrange's planetary 
equations to find the perturbation in the inclination of the orbit of an earth sat- 
ellite. Before doing this, however, let us make several simplifications: (a) omit 
the zero sum -ra 6 *+m 9 , (b) consider only terms for wMch j = 0 and £ = 2h, 
to be rid of dependence on the satellite eccentricity and mean anomaly, and (c) 
consider only terms for which q = 0 to be rid of small terms which depend on 
the eccentricity of the disturbing body. The total tidal potential now becomes 


-r =ti Cm* ( 2 -m) ! ^ /M 

^2mp0 '•‘impO (2 + m)! \ 117 \ / 


• F. (1*)F, ,(l)cos (2-2p)(M^ + cj*) + mf2*-mf2 + 5 , 

2mp ' ' 2ml 


mpO ' 
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Substitution of this potential into the Lagrange equation 
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and integration with respect to time yields 

^2mp0 


Ai, 


mpO 


VcMpU (I -e')''- sin 1 


(2-0,, ) m 

Om + ni)! 




cos |(2 - 2p)(w* + M*) + mn* - mJ2 + 


F (l'‘‘)F (1) 

[(2-2p)(aj* + M*) + mn="-mnj 


(2) 


We have assumed in the integration that the satellite elements, a, e, and I, as 

• • • 

well as the nodal rate are constants. We have also assumed that 

*, and I* remain constant. This last set of assumptions must be removed in 
the case of the moon when considering periods of time longer than a few years, 
due principally to the motion of the lunar node on the ecliptic. 

The tidal inclination perturbation is thus given by 

AI= ) AL . + > AIt n+AJ„ 

/ 2mp0 / ImpO 0 

ni,p iiifp 

Moon Sun 

/r • 

where A Iq is a constant of integration. 

We note in passing several features of eq. (1)* First, tides for which m — 0 
do not perturb the orbital inclination; thus we obtain no information on these tides. 
Second, the periods of the tides felt by the satellite are greatly lengthened 
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over their periods on the earth's surface, due to the replacement of 0 in the 
frequency by tl. Third, there is no way to separate the K ^ tide from the K,s 
tide ( P = 2110 ), or the tide from the tide ( 2^^ -2210), since in 

either case the periods are the same. A hypothesis is needed to do tMs, 


ANALYSIS OF THE DATA 


Sq. (2) may be conveniently written 
-Sp,„po I 1 Sin 


(2-2p)(co=’= 

(2-2p){co' 


+ M*) + - mJ2j 

+ M=-’) + 


for the analysis of the data. Here cos 

sin , and the empty brackets denote a factor common to both terms. The 

2mp0 ’ ^ 

coefficients c, „ and s. ,, are found from multiple linear regression, from 

which the effective Love numbers and phase angles may be easily computed. 

The regression analyses for the three satellites were carried out with the 
computer program described in Appendix 3. In each case the tides with the 
largest expected ampHtudes in eq. (2) were retained in the analysis and the 
smaller ones ignored. The theoretical amplitudes in the solid earth signal were 
conveniently used to pick out the major contributors. 

Table 2 gives the regression solutions. 'The effective Love numbers and 
phase angles for the K m , and K s tides are assumed to be the same; likewise 
for the and tides. The resulting incUnation curves and residuals are 

shown in Figures 1 - 6 . 
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Another computer run was made in which the effective phase angles of all 
but the Kj tide in the BE-C data and the Kj and S, tides in the GEOS-I and 
GEOS-II data (the major contributors in each signal) were set to zero. No major 
changes in the effective Love numbers (or in those effective phase angles solved 
for) took place. A similar computer run was made vdfh only the first 28 points 
in the BE-C data and the first 72 points of the GEOS-I data retained. Also, 20 
neighboring points in the GEOS-II data taken from the middle of July to the end 
of September in 1968 were removed and a run made. Again no major changes in 
the effective Love numbers and phase angles took place, indicating the stability 
of the solutions. However, the deleted points in the GEOS-H data were chosen 
not to disturb the basic signal, due to the spacing of the data. 

Examination of Table 2 reveals good agreement between the BE-C and 
GEOS-I results, but poor agreement between GEOS-H and the other two satellites. 
Also, there is good agreement between our determination of ocean tide param^ 
etprs and that of Eelsentreger, et al. (1976) for GEOS-I, but no agreement on 
the ocean tide parameters for GEOS-II (see below). The cause of this discrep- 
ancy is unknown; hence the GEOS-H results must be treated with caution. 

Most of the effective Love numbers given in Table 2 are smaller than the 
solid earth Love number k^ of 0.29 to 0.31. This is in general accord with the 
ocean tide models and charts given in Table 2 of Lambeck, et al. (1974). The 
only entry in that table which predicts an effective Love number greater than k, 
is the O^ tide chart of Ehetfich. (The S^ tide model of Bogdanov and Magarlk 


predicts an apparent increase in the effective Love number for GEOS-n due to a 

fourth degree term.) The phase angles are on the order of a few degrees, as 

expected, although some are negative and the BE-C find the GEOS-H 

effective phase angles are embarrassingly large. Also, the BE-C and GEOS-n 

S effective Love numbers are rather large and the GEOS-n P effective Love 
2 * 

number is rather small. These presumably reflect the inadequacies of the data. 
However, it will be noted that all three satellites indicate a remarkably small 
effective Love number for the P^ tide, apparently indicating a large ocean tide 
effect. If this effect is real it is highl;/ interesting, since the tide (the lunar 
counterpart of the Pj tide) shows no similar behavior. In fact, the effective 
Love number should be larger than according to Dietrich’s chart, as men- 
tioned above. Unfortunately, we have only one poorly determined effective Love 
number for the tide, and no mathematical models for either the or P^ 
tide appear to exist. 

SEPARATION OE OCEAN AND EARTH TIDES 

We now solve for the ocean tide parameters from the satellite data in the 
manner of Felsentreger, et al. (1976). To do so, we assiune that the tidal signal 
consists solely of the solid earth and ocean tides ; this lumps the small atmo- 
spheric tides (principally an S;, tide) with the ocean tides. 

Separation of the ocean from the solid earth tides requires knowledge of 
the solid earth tidal response. If we assume the solid earth Love number to be 
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frequency-independent, then from Appendix 2 we have 

in even 


in odd 

(3) 


in even 


m odd 


(4) 


®2mp0 “ ‘^2mp0 ^2nip0 *^2 2mp0 mp 


, = sin e ( 1 + kV) 


^(2mp0)2tn *^2i 


'^(2mp0)2m ^2in 


^2mp0 “ *^2nip0 ^2mp0 *^2 ‘=2mp0 ^mp^* *^2 ^ 


^(2mp0)2m ^2m 


^(2mp0)2m ^2m 


Where e « m<^», {Lambeck, et al. 1974), 0, being the solid earth lag angle 
and k * the second degree load deformation coefficient. 

The question arises as to what values to choose for k,, kj, and 
pick k^ - 0.3G and K' = ■^0.30, in accord with Longman (1966), and 0, * 0, in 
accord with a high Q { ~ 230) for the earth’s mantle (see, e.g. Stacey, 1969). 

The two rightmost columns of Table 2 give the resulting ocean tide amplitudes 
and phases. These are to be compared with Table 2 of Lambeck, et al, (1974) and 
Table 1 of Felsentreger, et al. (1976). 

The ocean tide amplitudes in our Table 2 have absorbed the fourth (and higher) 
degree terms in the ocean tidal potential, which Lambeck, et al. (1974) showed 
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to be important for some tides. The fourth degree terms could be solved for 
from observations of two or more satellites, since these terms cause an appar- 
ent variation in the second degree effective Love number from satellite to satel- 
lite (see Appendix 2). We have chosen not to solve for them here, since the 
differences between the effective Love numbers tend to be several times larger 
than those predicted by the ocean tide models and charts; aiid there are no clear 
trends . 

We may also make a rough estimate of k, from the ocean tides charts and 
models, in the manner of Lambeck, et al. (1974), For tMs we set 


d, +Z (1+k;) 


sin 

^(2mp0)2m ^2m 


^(2mp0)2ni ‘^2in 


m odti 


and substitute the ocean tide amplitudes and pirtses in Table 2 of Lambeck, 
et al. (1974) into the right side of this equation and our values into the left. A 
simple average of all charts, models and satellites yields k, « 0.27. Perform- 
ing the same calculation without the GEOS -II results gives k, =« 0.26. 

The effective Love numbers and phase angles of the wmll -determined tides 
permit an estimate of 0, as well as k,. Using the K, tide from BE-C and the 
and S, tides from GEOS-I and GEOS-n, in connection with the and 
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tides from Table 2 of Lambeck, et al. (1974), we find from eqs. (3) and (4) that 
K, w 0.29 and 0, «sl.5“. Once again an arithmetic average has been taken. 

From 02 determine a rough estimate of Q for the mantle, since (see, 
e.g., Stacey, 1969) 

Q W ! - as ! 

^ tan 6,200 -^2 

yielding a value of 19. Omitting the GEOS-n results yields k, « 0.28, 0, ^==^4.2°, 
and Q 7. 

Only a simple arithmetic average is taken in the computations above because 
we have treated the ocean tide models and charts as ’’observations” without know- 
ing how to weight them. Further, the quality of the satellite data prevents us 
from reading too much significance into the numbers (such as the extremely low 
Q values). Hence the values computed above indicate only the plausibility of the 
hypothesis that the ocean tides are responsible for the variations in the effec- 
tive Love numbers and phase angles. 

TIDAL FRICTION 

The evolution of the earth-^moon system is dominated by the second degree 
terms in the total tidal potential. In fact the secular rates of change of n^ , 

> ®M > I s ’ ^ written 

m pq 
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where X is one of the above-mentioned quantities and the right side of the equa- 
tion has the form (Appendix 1) 

t^^2mpq ~ * 2mpq®2nipq 

The values of „ are given in Table 5. Note that the s are determined 
directly from the data by the regression analysis, apart from possible corrup- 

'•n . 

tion by the fourth degree ocean tide terms. 

The [X] 2 ,„pjj may be conveniently arranged in a table, with the entries to be 
filled in as reliable data become available. Our owm tentative values are found 
in Table 3, These were computed from the weighted average over all three 
satellites of each Table 4 is similar to Table 3, except that the GEOS-H 

data have been omitted. 

Our value of n^^ » -29 ± 15 arc see/(100 yr)^, taken from Table 3, compares 
favorably with Lambeek’s (1975) -34 ± 5 are sec/(100 yr)2, computed from ocean 
tide models and charts. It also agrees well with the astronomically determined 
values, which range from -37,5 ± 5 to -52 ± 4 arc sec/(lG0 yr)" (Lambeck, 1975). 
The value of -17 ± 20 are see/(10G yr)“ from Table 4 is barely half the 
current estimates. 

Our value of E as -3.6 ± 1.8 x 10 *®erg/see from Table 3 is also in good agree- 
ment with current estimates. It is bracketed by Miller’s (1966) value of -1.7 x 
10 ergs/sec, determined from energy dissipation in shallow seas, and Lambeek’s 
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(1975) value of -5.7 x 10^^ ergs/see, determined from ocean tide charts and 
models. The number quoted in Table 4 is close to Miller’s value, which probably 
represents a lower limit on the energy dissipation. 

On the whole, it would appear that Table 3 agrees better with modern esti- 
mates than Table 4. However, Table 3 Indicates that the effective phase angle 
of the S 2 tide is negative, which would speed up the earth’s rotation. This is 
reminiscent of Holmberg’s (1952) hypothesis that the atmospheric tidal torque, 
which tends to speed up the earth, balances the frictional torque, which tends to 
slow the earth down, thus keeping the earth at a steady rotation rate. The data 
given in Table 3 of Lambeck (1975) does not indicate a negative effective phase 
angle for the tide, however. 

Our data on the tide may be used to estimate the time scale of the tidal 
evolution of the earth-moon system. Taking (see Appendix 1) 




'M 


^2200 “ 


3M,-a 


-.11/2 


^* 220^^M ” ^^2200 ^2200' 




multiplying each side by a^^^" and integrating with respect to time juelds 


ji3/2 _ ai3/2 ^ i3/6VG(M|, ) 


nn 


M, 


Ri [F 


220 " ^2200 ^2200^*^ ^0^ 


assuming the inclination. Love number, and phase angle remain constant. At the 

« 

early time t^ when the moon was close to the earth a^ « 0. Hence the tidal fric- 
tion age of the earth -moon system is approximately 
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T« 6/13 




V G(Mjj + ninj) m,^j 1^^220 "^^2200 ^2200 

0.0469 


^2200 ®2200 


X lO*-* 


yr. 


Tlie above equation gives T = 1.4 x 10^ yr using the weighted average s^qq 
value for all the satellites, and T = 2.4 x 10® yr using the BE-C and GEOS-I data 
only. These times are obviously much shorter than the age of 4.6 x 10® yr for 
the earth and moon determined by radioactive dating. A way out of Lliia Lixrre 
scale difficulty has been pointed out by Lambeck (1975), He notes that the 
tidal torques may have been much smaller in the distant past when the oceans 
did not flood the continental shelves, and that the ocean basin configuration was 
much different than it is now, due to continental drift. 


DISCUSSION 

One might wonder whether the variation in the effective Love numbers is in 
fact due to the ocean tides and not due to frequency-dependence of the solid earth 
Love number. Alterman, et al. (1959) show that k, does vary with frequency. 
However, the predicted change in Love number over the tidal frequencies is 
much smaller than that observed here. 

Another possibility is core-mantle resonance. The theories of Jeffreys and 
Vincente (1957) and Molodensky (1961) both predict a variation in Love number 
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near the diurnal frequency on the same order of magnitude as that observed 
here. These theories would not explain tlie observed differences at frequencies 
away from the diurnal frequency; but core-mantle resonance cannot be ruled 
out as a possible contributor to the Love number variation. 

Our results lend support to the ocean tide hypothesis of Lambeck, et al. 
(1974), and indicate what can be done when more data become available. They 
also indicate directions for further study. Regularly spaced and densely packed 
laser observations, for instance, would give good information on the short- 
period tide. Observations of a low-inclination satellite would also increase 
the amplitude of the signal, thus determining it more reliably. If the effec- 
tive Love number of the Pj tide is really as small as our data suggests, then a 
mathematical model of the Pj ocean tide would be of great interest. 

NOTATION 

a*I=*e*M+f2?*=n; * orbital elements of disturbing body (sun or moon) 


aleMfico orbital elements of satellite 


C 


moment of intertia of the earth 



ocean amplitude 
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1 1 1 1 ! 

effective tidal Love number 

1 

G 

universal constant of gravitation 


's 

solid earth Love number, degree 2 

* 

k; 

load deformation coefficient, degree s 


m* 

mass of disturbing body (sun or moon) 


M,. 

mass of the earth 


27rtifl’ 

tidal argument 



radius of the earth 


E 



^Cmpq 

effective tidal phase angle 

> 

^£mpq 

solid earth tidal phase angle 


+ 

^sT 

ocean tide phase angle 
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Q 


rotational speed of the earth 


pj, average density of the earth 

density of sea water 

0 , solid earth tidal lag angle 
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APPENDIX 1 


We discuss here th.e rates of charge of the lunar orbital elements and the 
orientation and spin of the earth under the influence of tidal friction, and the dis 
sipation of tidal energy. 

We treat the earth -moon-sun system as two separate systems, namely 
earth-moon and earth-sun. Strictly spealcing, this is imp ermis sable, since we 
are dealing with a tlaree body problem. However, the errors made in doing this 
are small compared to the uncertainties in determining the effective tidal Love 
numbers and phase angles. 

In the following equations we follow the formalism of Kaula (19G4), using 
the sign convention of Lambeck, et al. (1974) and the corrections for the two 
body problem of Rubincam (1975, Appendix 1). 

The secular rate of change in time of the lunar semimajor axis is given 
by (Kaula, 1964) 

6,m,p,q 


where 







r22+1 

,,2Jd+3/2 


(2 - ni)! 
(2 + m)! 


(2-5o,„){2-2i3 + q) 


^ ^'empq S'” 
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The moon's mean motion is related to by 


VG(Mg + m^,) 

'\m 3/2 


so that 
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dt 
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The components of the rate of change of the lunar orbital eccentricity are 


given by (Kaula, 1964) 
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Let us call the angle between the plane of the lunar orbit and the ecliptic J. 
Then we have approximately (Rubincam, 1975, Appendix 1) 


Rf'-' (C-m)! 

[jt(, = \/G(M, + m., ) — 

'Cmpq I.-. M ,, 2C+5/2 


M,, a 


M 


(y + m)! 




(C-2p)cos 1^, -m ^ 

, ,, f^Cpq ^ ^Giiip ^ ^ ^ '■^Cinp<i ^Gmpq ■ 

(1 -e-,)"sin I_^,, 


We now turn our attention to the earth. Note that both lunar and solar tidal 
friction disturb the earth. This is in contrast to the treatment above, where only 
lunar tidal friction disturbs the moon’s orbit. The tidal bulge raised by the moon 
is geared to the lunar motion, but the solar tidal bulge is not. 

From considerations of the conservation of angular momentum the rate of 


change of the earth's equatorial tilt to the ecliptic is given by (Rubincam, 1975, 


Appendix 1) 



Gmpq 


R,-C0 \a*/ (C + m)! 




(G - 2p) ~ m cos I’*’ 
sin i* 




and the slowdown in the earth's rotation by 
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Let us now derive the rate of dissipation of tidal energy in tlie earth. 


The orhital energy in the center of mass frame is 


E 


0 




and the rotational energy of the earth is 


E^ = ViC0- 


Therefore the total meehanical energy is simply 


E= Eq + E^-- 


GM,.m=* 


2a=" 


+ ViCO-. 


Thus the rate of dissipation of energy is 


E — 


GM,..m* 


^ +C0 
dt 


di' 


24 


I 


From the previous equations we derive the relatively compact expression 




(g - m)! 

(.e+m)! 


- Son. ) 


• [(2 - 2,p + q)n* - mO ] 


It will be noted that the above equations all have the form 


gnlpq ^empq ®Smpq 


where X is the quantity of interest, sin and Pg,„pt,is a multi- 

plicative factor whose value is easily computed. Table 5 gives the Pgp,p^ for 
2=2 (the degree which dominates the evolution of the earth -moon system) and 
q = 0. In the computations for p 2 ,„po lunar inclination functions Fg^^p (I 
which change with time due to the motion of the lunar mode, axe assumed to take 
on their average values. Also, GjppCe*) is taken to be 1. 

Table 5 gives us some insight in determining which tides are important in 
changing what quantities. It is clear that the Mt tide governs the rates of change 
of n a„, 0, and E. The O „ K^’, and M, tides appear to be about equally im- 
portant in changing J and I^. The table gives the appearance of the M 2 tide gov- 
erning the change in ; but it is probably the N 2 < 2mpci = 2201) and L, ( £mpq * 
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220{-l) ) tides which dominate the change in the lunar orbital eccentricity. These 
tides are not given in the table. 
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APPENDIX 2 


We wish to show here how the ocean tides combine with the solid earth tides 
to give the total tidal potential, as explained by Lambeck, et al. (1974). 

The potential of a particular ocean tide constituent n is given by (Lambeck, 
et al., 1974) 
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n 


4irCR^. 
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L-cos. 


(s - 2u)co + (s - 2u + v)M + t(n - 0) ± 27riifT + e[^ 


s - i even 
s - 1 odd 


Note that a spherical harmonic in the tide-raising potential excites many spheri- 


cal harmonics in the ocean tide response; hence the summations in the equation 


above. 


We may simplify^ the expression above by getting rid of small terms which 
depend upon the satellite’s mean anomaly M and eccentricity e by setting v = 0, 
nn d s ^ 2u + V = 0 (which implies s is even and u = s/2). The equa- 
tion then becomes 
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Let us now convert the tidal frequency into orbital elements with the cor- 
respondence 


-27rnfT = (£ - 2p)co* + (£ - 2p + q) M* + ni - fJ) - niTr. 


The correspondence is not exact (Lambeck, et al., 1973) and more detailed con- 
siderations must be given if the period of time under consideration is more than 
a few years. We now have 
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s - 1 even 


s - 1 odd 


where the subscript n has been replaced by ( Cmpq). 

Let us now get rid of terms which depend upon the earth’s rotation 0. The 
only way to do this is to pick the plus (+) in S and set t = m. Then 
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Also (47 t/ 3) ppRp* so that we may write 
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where the empty curly brackets mean 


= (2 - 2p)co* + (2 - 2p + q)M* + mfl* - mf2. 


Compare this to the solid earth tidal potential, which for 2 = 2; m > 0; 


q = 0 and small eccentricities is 
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Assuming that the total tidal potential is composed of the solid earth and 


ocean tides only, one may show from the expressions for the ocean and solid 
earth tidal potentials that 

‘^2mp0 ^2mp0 ” ^ ^2mp0 ''' ^mp *^2^ 


'^(2nip0)2m^‘” ‘'2in 


+ \ 111 even 
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(2 + m)! /M,A Ip \/a*V > 
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(2-m)! \m*l\p^J\Rj F2,„p(l*)R^ 


The term in brackets at the end of each of the above equations for d^ „cos 6^ 

^ 2mp0 2m 

and d,^pp sin represents the amount by which the second degree effective 


Love numbers and phase angles will be corrupted by the fourth degree ocean tides , 
if the fourth degree terms are not solved for. 

Care must be taken when solving for the effective Love numbers and phase 
angles of the K‘|' and K® tides. The reason is that the and e'^'are given for the 
Kj = K‘1’ + K® tide, and not for the K^'and K® tides individually. Some hypothesis 
is needed to decompose the Kj into the K*'* and K® tides. Similar statements hold 
for the and tides. 
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APPENDIX 3 


A, PURPOSE 

Tlie computer program solves for the effective frequency -dependent tidal 
Love numbers and phase angles for the whole earth (solid earth + oceans + 
atmosphere). It does it by performing a multiple regression analysis on the 
observed tidal perturbation in the inclination of the orbit of an earth satellite. 

It also solves for ocean tide amplitudes and phases from the Love numbers and 
lag angles found from the regression analysis, given an assumed tidal response 
for the solid earth. 

All tides are considered to be of second degree only, and only tides which 
depend on the first power of the orbital eccentricity of the disturbing body are 
retained; hence the subscripts on the inclination and eccentricity functions run 
as follows: g = 2; ro = 1, 2;p = 0, 1, 2; and q = -1, 0, 1, Within these limits, any 
combination of tidal constituents may be chosen to fit the data, giving their ef- 
fective Love numbers and effective phase angles. Also, the phase angle of any 
constituent solved for can be set to zero, resulcing in a solution for the Love 
number only for that constituent. There is no restriction on the number of dif- 
ferent regression analyses that can be carried out on the data during a run of 
the program. For example, the sample data cards given in the listing of the 
main program give two regression analyses. 
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In practice only the few constituents which give the largest signal in the 
inclination are fit. If these are not known, then a call to subroutine THEORY 
gives the theoretical amplitudes and periods of all the constituents in the solid 
earth signal, given an assumed Love number for the solid earth, such as 0.30. 
These can be used as a guide for picking out the important constituents. The 
subroutine also gives a sample plot of inclination vs. time, if this is desired. 

The precession of the lunar orbit around the ecliptic requires a certain 
amount of averaging to be done on the input data for the moon; this is explained 
in section D. Suffice it to say here that the motion of the lunar xiode and the 
moon in its orbit with respect to the earth's equator are average rates over the 
data span; the inclination of the lunar orbit to the earth's equator is assumed to 
be constant and equal to its average over the data span. If the data span is longer 
than four or five years, then a different program taking the precession of the 
orbit into account more exactly will be needed. 

B. THE EXPRESSION FOR THE TIDAL INCLINATION PERTURBATION 
The expression for the tidal inclination perturbation is 




^^2mpq ^^2nipq 


(2-m)! 
(2 + ni)! 


Gm* 
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cos 



(2 - 2p) (cj* + M*) + qM=^ + mr2* - mO + 6 
1(2 - 2p) (w* + M*) + c|M* + - mfi ] 


The program names of the variables in the above equation and in wliich sub- 
routine they are computed are given below. 


Variable 

Program Name 

Subroutine 

ci-. 

2mpq 

YKM(LM,LP,LQ),YKS(LMXP,LO) 

REGRESXOVNUM 

^2iii|jq 

EPS LNM( LM ,LP, LQ),EPS LNS( LM ,LP,LQ) 

REGRHSXOVl>JUM 

2 (2 - m)! 
(2 + m)! 

Bl(LM) 

BLM 

m 

YLM 

LOVNUM 

Gm* r| 

VgM,, 

CFNTM, CENTS 

COEFF 

I* 

XIMOON.X1SUN 

READMS 


AI(LM,LP),B{LM,LP) 

INCL 

1 

XI 

RDSAT 

P2mp 

C(LM,LP) 

INCL 

e* 

EMOON,ESUN 

ECCFUN 

2pq ^ '' 

GLPQM(LPXQ),GLPQS(LP,LQ) 

ECCFUN 

sin Iq 

SS 

MAIN 

CO* + M* 

XNDM.XNDS 

READMS,MAIN 

M* 

XMDOTNhXMDOTS 

READiMS,MAIN 


34 


I 


Variable 

Program Name 

SiibrouUiie 

h* 

OMGDOT 

READMS 

O 

Q 

RDSAT 

(2 - 2p) (ojf + 

+ qM* + mii* - mO 

ARGDTM(LM,LP,LQ),ARGDTS(LM,LP,LQ) 

ARGDOT 

w* + M* 

XNDM*(TT-TMS)+DELMN, 

XNDS*(TT-TMS)+DELSUN 

READMS, ARC 

M* 

XMDOTM*(TT-TMS)+XMEANM, 

XMDOTS*(TT-TMS)+XMBANS 

READMS,ARG 

n* 

OMGDOT*(TT-TMS)+OMEGAM, 

OMEGAS 

READMS, ARG 

n 

1 1 

Q*{TT-TZERO)+XNODE 

RDSAT, ARG 

cos 

ARGUM(LM.LP.LQ).ARGUS(LM,LP,LQ) 

ARC 

1 ) 

AIq 

ANS(1),XJ0 

REGRES,LOVNUM 


The regression coefficients are found by writing 


AL = d, cos 5, [ 

2mpq 2mpq 2mpq 


-d 


2mpq sin [ 


] cos 


(2 - 2p) (oj* + M*) + 


mi2 


1 sin (2 - 2p) { lj *' M*) + qM* + “ mJ2 



so that tbe data are of the form 

^ a.x,. + b 
j 

where the regression coefficients Uj are d^^^pq cos 52,„pq and -djjppq sin Sjmpq’ 
b is AIq, d 2 ,^pq and 5,^,pq can obviously be solved for from the regression coefficients. 
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C. BRIEF SUMMARY OF OPERATION 


The program begins by reading in the satellite data; the sun and moon data; 
TSTART, TEND and DT; the assumed solid earth Love number and load de- 
formation coefficient; and the solid earth tidal lag angle. Next, parts of the 
tidal inclination perturbation are computed. Subroutines TIMEl and CORRES 
are called, which are used in the plotting and printing of the inclination pertur- 
bation and residuals against time. 

The program then reads in JTHORY and JPLOT; if JTHORY = 0, then the 
regression analysis is done wihtout calling subroutine THEORY. If JTHORY = 1, 
then subroutine THEORY is called but no regression analysis is done. JTHORY 
= 2 does both. If subroutine THEORY is called and a sample plot of the solid 
earth signal is desired, then JPLOT should be 1. If no plot of the solid earth 
signal is desired, then set JPLOT = 0. 

The regression analysis begins by calling subroutine LASDTA, which reads 
in the observed inclination. Subroutine REGRES is then called, which reads data 
cards to determine what constituents are to be solved for. The comments given 
in REGRES tell how to pick the constituents. Subroutine REGRES computes the 
regression coefficients and their errors, prints them and calls subroutine 
LOVNUM. 

Subroutine LOVNUM computes the effective frequency-dependent Love 
numbers and phase angles and their errors from the regression coefficients and 
their errors. Then it computes the ocean tide amplitudes and phases from these 
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by subtracting out the solid earth tide. All of these results are then printed and 
the inclination and residuals vs. time are printed and plotted. 


D. SAMPLE INPUT 

The sample input is listed in the main program. Each data card (number) 
and each number on the data card (letter) are discussed below. 

1. The satellite name and when the data were taken are given. 

2. a. The semimajor axis of the orbit is 8.07291 x lO® cm. 

b. The eccentricity of the orbit is 0.G726. 

c. The rate at which the node progresses along the earth’s equator is 
-2.246536 degrees/day. This figure is obtained by noting that on Feb. 7, 
1966 near the beginning of the data span the node position was -109. 082 02 
degrees, and on Oct. 26, 1967 near the end of the data span the node posi- 
tion was -75.41322 degrees. (All times are 0 hours UT.) Hence the 
rate is 

-1440 -t- 109.08202 - 75.41322 ^ _2.246536 degrees/day 
626 

where the -1440 = -4(360) indicates four complete revolutions of the 
node. 

d. The inclination of the orbit is 59,38053 degrees with respect to the 
earth’s equator. 
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e. TZERO is Feb. 7, 1966 12 hrs UT, so TZERO = 38.5. 

f. The node position of the satellite at time TZERO was -109.08202 de- 


grees. 

The American Ephemeris and Nautical Almanac (Nautical Almanac for short) 
and Figure 7 are used for obtaining the numbers on the next two data cards. 

3. a. The positions of the moon and sun are given for Feb. 4, 1966 in the 
Nautical Almanac, near the beginning of the data span; hence TMS is 
taken as TMS = 35.0. On this day the position of the lunar node on the 
ecliptic was = 60.8500 degrees, according to the Almanac. The posi- 
tion of the lunar node on the earth's equator is found from the formula 


tan 


sin J sin Q 

sin \ cos J + cos I sin J cos SI 


where J = 5.1453964 degrees = inclination of the lunar orbit to the 
ecliptic, and I = 23.44319 degrees = inclination of the ecliptic to the 
earth's equator. Solving the equation gives SV'"= 10.1770 degrees, 
b. DELTA! is M* (measured from the earth's equator); so 

DELTAl = = {( - n + c 


where (( and SI are given by the Almanac and c is the arc sho^vn in 
Figure 7. c can be found from the formula 

cos 1 ^ cos I* cos J 

cos c = ^ — .J, , 

sin I* sm J 
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where I*, the inclination of the lunar orbit to the earth’s equator is given 


by 


cos [* = cos 1 cos J — sin I sin J cos O 


We find that on Feb. 4, 1966 I* = 26.3140 degrees, c = 51.6092 degrees, 
(( = 115.2340 degrees, and 60.8500 degrees. Hence DELTAl = 

105.9932 degrees. 

G. The nodal position of the sun is always 0.0 degrees. 

d. DELTA2 = cc* + M* for the sun is 313.7163 degrees, according to the 

Almanac. 

e. The speed at wMch the moon moves around in its orbit averaged over 
the data span is XNDMl = co* + M* = 13.183848 degrees/ day. This is 
found by computing the lunar position on Oct. 27 , 1967 at the end of the 
data span. We find that SI* = 5.0417 degrees, I* = 28.1025 degrees, 

c = 22.9427 degrees, and cu<= + = 131,8174 degrees. Subtracting 

this last figure from cu* + for Feb. 4, 1966 and making allowance 
for 23 revolutions of the moon gives 

XNDMl = 13 1.8174 - 105.9932 ^ 13^3^33348 degrees/day. 


f. The average speed of the moon’s node on the equator is 


OMGDTl = 


5.0417 - 10.1770 
630 


-0.00815 degrees/day. 
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4. 


a. The inclination of the moon’s orbit to the earth’s equator averaged over 


the data span is 


XIDl = 


26.3140 + 28.1025 
2 


= 27.2082 degrees 


b. The inclination of the solar orbit to the earth's equator is 23.44319 de- 
grees . 

c. TMS is 35.0 days (Feb. 4, 1966). 

d. The position of the lunar mean anomaly on Feb. 4, 1966 is (( - r ’ = 
331.5588 degrees, according to the Almanac. 

e. The solar mean anomaly on Feb. 4, 1966 is 31.3590 degrees, according 
to the Almanac. 

5. a. TSTART is 37.0 days (Feb. 6, 1966). 

b. TEND is 666.0 days (Oct. 28, 1967). 

c. DT is 2.0 days. 

6. a. The solid earth Love number is taken to be 0.30. 

b. The load deformation coefficient of the solid earth is taken to be -0.30. 

7. The solid earth tidal lag angle is taken to be 0.0. 

8. a. Both the theoretical amplitudes and regression analysis are wanted, so 

JTHORY=2. 

b. A plot of the theoretical tidal signal is wanted, so JPLOT = 1. 

9. a. There are 142 data points. 

b. TDAY is 38.5, since this is the date of the first data point and the time 

on the following cards is measured from the first data point. 

40 



10. a. 

b. 

11. a. 

b. 

12-151. 

152. a. 

b. 

c. 

d. 

153. a. 
b- 

154. a. 
b- 

d- 


The tidal inclination perturbation of the first data point is -1.1594475 
arc seconds. 

The time of the first data point is 0.0 days. 

The tidal inclination perturbation of the second data point is -0.89708115 
arc seconds. 

The time of the second data point is 13.0 days after the first. 

, The rest of the data cards. 

The name of the satellite is GEOS I. 

There are 142 data points. 

Tliere are 68 independent variables + 1 dependent variable = 69 vari- 
ables. 

Two subset selection cards follow. 

Eight variables are chosen for this regression analysis. 

1. Fit the Kj tide (9 and 10), the K,tide (27 and 28), the Pj tide (39 and 40). 
and the S^ tide (55 and 56) for effective Love numbers and phase angles. 
Seven variables are chosen for this regression analysis, 

•c. Fit the K^and s; tides for effective Love numbers and phase angles 
(9 and 10, 55 and 56). 

■g. Fit the M 2 , K,, and tides for effective Love numbers but no 
phase angles (which are forced to be zero). 
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COMPUTER PROGRAM LISTING 
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main PRU3»AM 


C THI5 JHLVt: Fnn TMt£ cPPHCtlV£ FR£0' tENCt - CEPtN.'.'fcNT TIDAL LOVE 

C NUMOcP^ AND UA<3 AM3L=$ FOP the *1HaLE EA«TM flY CARRY.IKO A MULTIPLE 

C Rl3RES*IDN analysis on ThC. tidal PERTUPEATICN IN THE IKCLJNATICN OF Tt-£ CPUIT 
C OF AN EART H SA'^FLLfYs* ji AlSC ?0LV“3 FDP OCEAN TjnE APPLlTUOfcS AND Ph *?6S 
C PNOM TMs lP'/c. NUA^lEhS anC LaG AtJOLES FCOND« OlVEN Tt*r TIDAL RESPCNEE OF T>f 
C :aULlD EARTH. 

C 

C ONLY i-ND uEGREI TIDES APL CONSIDERED* Aii WELL AS CNLY TERMS WHTCt* CEPENO 

C ON THL rlf;>T POWER OF ORBITAL rCCLNTPIClTY HENCE THf. SUESCRtPTS ON THE 
C iNCulHATION AND cCCtNTHlClTV FUNCTIONS RUN AS FOLLOWS L = if* M= 1 , 2 , P=0*l,2 . 

L AND 0=-lfC*l • (SINCE THE COMF(J'S» C AUNOT HANDLE 2£PC SUBSCfilPTE* IN THIS 
C PRO-SWAM THEY PIJ^^ As F'lLLC.VS LM=M, LP=“-Hl, AND LC=CH2 .) .^ITHlN THESE 
"C^'LTM'fTs* "Any CCMrt INAY Jhn' of yicaL CCNETlTUSNrs MAY H£ CHC55N FOP FITTING TO 
C The data- 2m Pi^ACYIce* only THP few constituents which GIVP the largest signal 
C IN TH=. INCLaNATICN ARE FIT. (THE XMPCRTAfJT SIGNALS MAY CE FCUNO BY C-ALLING 
C siUrJMOUTI NE TH=:CPV.» 



— ^ — 
c 
c 
c 
c 
c 
c 
c 

DkIEF Dc.O( 

ahG • < 

ARGOOT - < 

At<GLST “ ( 



“OL7-r - < 


c 

COEFF - 1 


c 

( 


c 

CUP R Fa - 1 


c 

1 

— ■ ■ 

t 

OATA ~ ( 


c 

sCCFUN - ( 


c 

INCL - 1 


c 

iiAur — • 


c 

LAG.- 


c 

LA3DTA - 1 


c 

LOVE 


c 



c 

LOVNUH - ' 


c 

1 


c 

plutcp - 1 


IN this program “fir^AL INCLINA72CN PEPTyqfATION 


□kILF Oc.3CmPTlLN 'JF EACH SUtlPCJTlNE 

AHG • COMCLTES the COSINE. Tgp^/ the expression FCR ThF tIP 
ARGOOT - CUMPHTE5 THE OFPIVATIVE OF THf ARGUMENT OF THE COSINE TERM IN 
THE TIP 

AkGLST “ CCMPU^Et SINES AND CHSlNSS TO BE UScC JN SLIrKCUTTNE DATA 

ou-r CD’iPLTtS The FaCTOr I al" TL”M In THE'T'IP 

COEFF - CCMRijTES THE CUEFFI c 1 £nt i*,j TME TIP CF THE FCPM 
CCNoTANT/( A** 7/2 J 

COPREa - UEEC only RY FIJOROUTINE olUTLF TC OFTEPMJnc ^HEThEF AN 

EXPERIMENTAL DATA nQINT SHIULC at PLCTTEU AT A daPTICULAC TImE 
data - CCMPUTE5 VALUES TO «?£ L2r.O IN REGRESSION ANALYSIS 
sCCFUN - CCMPUTES ECCENtpICI^Y FI-nCTIONE 

INC L - CCMFLTZS INCLINATION F'JNCTICNG _ 

‘CKZT' - '^EAD'E' IN SoLiD 'EAPTH LAC ANGLE 

LAG.- - THE I'h'EO'.JENCY DEPENDENT LAG ANGLES INITIALLY ECUAL TC 2EFO 

LASDTjA - RtADE Ifi THE CBSERVED TIP 

LOV£ - tpaDS in Th£ ‘iCLir> EAR'^H LOVE NLfMGEP AND LCAC DEFGPMaTSCN 
CLcFFTCIFnT 

LUVNUH - TAKES THE OUTPUT OF PEGPcS ANC COMPUTES FRtCUENCY-DtPENDEM LCVF 
WUMUERj* LAG ANGLE?* AND OCEAN TIDE PARAMcTER?- 
PLUTCP - PLTlb CUT TH? JIP« EXPF FI vJ-nTAL OaTA PCINTS ANC THE RcSICUAUE 

ON THE COMPUTER PApCP ‘ ' 

HD3AT - PEACE IN SaTcLLITF ELEMENT^ 

KEADML - READL iN SUN AND MOCN PCSITIDNS* FATES 
kEGRHS; - PuMFCTms RFGRCSSinN ANALYSIS 

7HEJRY - CDHpuVCb. AMPLITUPFC* PCHICDC. rjp tHE 5CL IP EARTH TICAL SIGNAL 
-TIHEJ - CCMPUTcE TIME VALUES FOO ARRAY TIL) 


notation 


satellite 


Q 

oATI J) 
Sri 

T7EPQ 


SEMIMAJOr AXIS IN 10»*3 CM 
iNCLlNATlaM FUNCTIONS 

CRBITal ECC^NTfilCItY 

RATE OF change of NODE POSITION 

ALPHANUMERIC INFOS. VAT ION * SUCH AS NAM£ OF SATELLITE 
SINE OF SATELLITE INCLINATION 

TIME (MrAFUREO FROM OEGINNING OF The VcAR) whSN THE 

node hao Value xnode 
satellite iNCLf naticn 

PCoITION DF NODE AT TIME TZERC 


MOON AND SUN 


lOOB QUA^m 
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1 


1 


I 


C ACGUM(LH»LP« tQ) 

C AMGU5(LM«(.P»LQt 

C 


C AROUTMiLHfLPfLC) 



C AtCLMfLPI 

C 

C 

C CFnTm 

t CFNTS 

C 

c oeuMK 

c ' 

C OfciLiUN 

c 

C cMCON 

C ESUN 

C 

C EPSi.NM(L4*LP *LO) 

C cPiLNS<LM,LP .UQJ 

c 

C OLPOMCLP*UOI 

C GLPQSU.P«4.0) 

C 

C OMCGAM 

C OMFOAt 

C 

C _ XIMCJK 

C xlS'JN " 

c 

c XMOOTk 

C XMOOTS 

c 

C XMOONtLM.UPtLO) 

C XSUN (LMiLP *L C) 

c 

X-40M 

c 

C XNJik 

C 

c VK«? LM#LP*LQ) 

c 

c VI*.:: < LP«LP#LO ) 

c 

^ 

C 

C UATA 


CC61NE ▼EOM IN THE TJP 

FATc OF CHANGE OF ARGUMENT OF CCSIKE IN THE TIP 

INCLINATION FUNCTIONS OF MCCN 

TERM IN TIP OF ^OPM CONST ANT /( #• • 7/2 I 

IA»G'JMEN.T QF_PERIGSii ♦ MEAN ANCMAm.._Cf- MCDh AT TIME.. 
TMS* MEAEUREO MITH RESPECT TO EAPTh*5 lOUATOR 

feCCENTRtCiTT OF LUNAR ORBIT 

EFFcCriVE FRFOUENCT-OEPENOENT LAG ANGLE FOR LUNAR TIDES 

ECCtNTRICI'^Y FUNCTION PCP MOON 

PCSITICN OF LUNAR NODE ON EAPTH*S cCLA^OP 

inclination of lunar ORBIT TO t*«Th*S iQUATCR 

RATE OF change of LUNAR MEAN ANCMALV 

LUNAR PART OF TIP 

PATE OF CHANGE CF (ARGUMENT OF PEtlCtE ♦ MEAN ANCMALVI 
OF MOON, ^ITH RESPECT TO THE EAHTH'S FOUATOP 

tFFtCTIVF FREQ lENCV-OEPFNCENT LOVE NUMPEP FCfi LUNAR 
T IOF«f 

tlM^'OF'MOdN AND aUN P03ITlbN!T 



Mc<P - NUMBER OF CA'^A POINTS 

TUAV - CORHtCTlON YO TgXP(J> TC C IVE iT iN CAYS FPCM THE 

HEOlNNlNG OF THE YEAR 

TFXP(J) - YIM£ OF OnSEPVATtnN XIEXP(J) 

XlExPrXJ" - CBFEPVEO TIP AT TIMt TExP(j) 

HiSCcLLANEOUS 


- INTERVAL OETXFEN SUCCESSIVE VALUE;? OF T(L» 

- total N1.IMOEP OF VALUES IN ARRAY T(L) 

- TIME VALUES BE^aFEN ANO INCLUDING TSTAOT AND TENf> 

- ENDING time OF TATA SPAN 

- STARTING TIME QFU^aTA span 

- CONVERSION faCYCR fROM -jEGREES TC iJADIANS 

- FACTORIAL PARY qF THE TIP 

- SOLID earth LCvE number (SECCNO CE5RFE) 

- LOAD OEFfRHATICN COEFFICIENT (SECLNC DEGREE) 

- AOOITIVF CCNF.Tan^ in TIP 

- LAG ANGLE CF SOLID EARTH TIDE 

- THt COMPUTED TIP AT TINE T(L) 

■- IN^EX ahicm <FCPS TCACic OF T I VE In SupPCuTInE Ca^a 

- OETcJJMINLf. AHETH'R a data PClN*^ bHCULC tic PLOTTED AT 
TIME T(L) 

- CCUTAIN" alphanumeric NAMES CF TICES (M*M2, ETC.) 

- USED IN SUPPOUTfNE LOVNUM AS A CCWAECTICN TO The 
LUNAR TIOFS TO GET THE PFtPER AMFlITUDE OF TH£ 

LUNISOLAC TIJEO 


c 

c 

c 

C 
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c 

L 

sAMPLt iNPUT (CU.UMK 1 nF 

THE XNPJT 

starts XN 

CTLUWK 3 HfcRE) 


C 

C 


( SAT( 31 

* 3-1*14) 



( format 

X3A6.A21 

c 


GEOt * t 

> 9€*-l 967 

DA*A 




c 

4 

c 

01 

Xll 

T2CWO 

XNOOEl (FCCmaT 

8FI0.5 1 

c 

d.07^91 

a* 0726 

-2*24c£3b 

59.38063 

38.5 

-1 C9.C620 


c 

c 

CME :*A 1 

delta 1 

OMcGAif 

TELT AT 

XNJWl 

CWGCTl IFJWWaT 

8F1 C.6 ) 

c 

13.1770 

105*9032 

0 * 0 

J! 3.7X63 

13.18385 

-0 .0C9I 5 


c 

c 

XlDl 

X1C2 

TVS 

XMC ANl 

XME AN2 

1 FORMAT 

9Fl0t5) 

C" 

■2T.2082 

2J* AA319 

35.0 

331 .EfFA 

31 .3590 



c 

c 

TSTART 

TtND 

OT 



( fllwat 

3F10.51 

c 

J7 *0 

660* C 

2 * 3 





c 

c 

XK 

XKP 




( FORMAT 

SFIO.S) 

c 

c 

0 * 10 

-0.30 







XLA^ 





( FC RM AT 

6F1 0 .5 1 

c 

3*0 







c 

c 

JTHUN Y 

JPLOT 




C FORMAT 

2X5) 

c 

2 

1 






c 

c 

c 

ME KP 

TO AY 
3S.S 




( FORMAT 

X5.F1 C.5 I 

r 

c 

xlEXP(J) "txP<J) 



( format 

2014.81 


OOO J 

0073" 


J3C. 13 

cTC. 


D 0? 


po ,PW J 

uFOj 1 


N M NS. 

*? 


N (:^AVF<J)« 
dk 9l02TibJ94055*t 
7 7 1021 


(F-t<MAT 36121 


T5fT 

:«N 
I 3N 
ISN 
I&^ 
1 SN 
I £N 
-T5K 


I SN oo: I 

ISN 0012 
ISN OCl 3 
ISN 0 : 1 A 
ISN 0 01 t 


0004 
000 = 
0006 
0007 
0006 
C 0C4 


I SN 
I SN 

» : N 
1 SM 
ISN 
1 *iN 

^rTTv 
I £N 
! SN 
I SN 
1 iN 
ISN 
ItK 
I SN 

TTN 
I 5N 
ISN 


•: 31 * 

031 7 

ooie 
c 3: 9 
002 0 
0 32 I 

-0032 
002 3 

302a 

C 32 e 
0 02 c 

C02 7 
0 02® 
^ ors 

03TC 
C 031 
00 32 


4 MPUICIT PF AL»4 (A-»-*0-2) ^ ^ ^ 

DT^rNnCN rPSLNM 12 »3tJ)*En^l.N£(2#lf )>• APOUMI 2 • 3 • J 1 ■ APt,C- ( 2t 3* - ) • 

1 APGCTMI i, 3, 3 ) • APGOTSI 2# 3* 3» 

JlMtNilCN hi (2J . XMCQM 2*3.3 1 • XSUM 2 .3 . 3 I . YK^ ( 2 .3 . 3 ) •YX^( 2* 3* 3) 
lIMLNS ICN TF XP <;'00 » * *IEXPI 200 ) *T{700 ) *Xt NCL ( 730> ACKI 00 1 

OlMtN* ION C< 2*3l*A*l(2#3)*fi(2.3) t^LPOMI 3*31*GLPQo !•-> 

COMMON/SCX A/ EP3LN«^ * cPSLNS * APGU*' • AOGUS « APOOTP .AOCuifS 
CJMMON/bLKfl/ ei*>«X*Cf'KTM,CFNTS* SU*** XMOON *XSJX*YK*»*YKS*XKP 

CCjMMCJN/r>»_KC/TEXP*XlPxP»T* kINCL.NTPACK . . wT 

COVVrTN/PLKO/ XlAO.O .T WE GAM* 2rLWK, CMr»5 xfi, OEL sUN *5S*XIC*A*t#xX*XSC0t* 
1 T 2_Ru. XI MC'JN.XI SUN.TMS* XNDM * X NO S • ^^ * XWt ANM, X WE AN 5 • X WCCT M, X MOOTS* 
d. CMGCCT*MCKP 

C JM*^CN/f3LXG/C *Al •t3*GLPQM*GLPO'“*»KOtX 
r = 3* 14 15 J2 S3 AD COO 
XNOS=< 0*M‘iS6A7OC)»F 
XMDQTM = ( I l*0e499200» *F 

XMDOTS *( F . 

PfAD it; ’TT^ r DATA* CJ3VOUTF PAPT > THF TIP 
CALL FDiAT 
CALL PsAOMS 

RcAJ (5.1) TSTART.TEND.OT 
FOP MAT ( 10 *0) 

A«ira (f-.Yi 1 syapt,t€no*ot 

Fi>P‘(Al (////////•’ 3X.*TsTAP^s**Pl)#>*5X*»TENC a* .Fl^.L. .X. CT* • 

4 Flw.‘ .&X. •(DAYS)* ) 

CALL XOVL 
CALL LACl 
call LA -.2 

j;,sOS!N( x: ) 

CALL eCCFDN 
CALL INCLC Xi *C) 

CALL INCL( X I MCCN* A1 ) 

CALL 1NCL(XI 3LN*0) 

CALL nuM 
CALL AS^ODUT 
CALL CCcrF 
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I IN 

OC • 

CALL 1 IM ifTSTAMT.TtND.DT.MTI 

1 'N 

0 0 54 

(NOLXsO 

1 5N 

00 iS 

MLXPsO 

I S-N 

0 Ole 

call cuphcsi ts^ apt, tend* OT .M t XP ) 

1 SN 

0 C 3 7 

RLA7 (5,2) JThO*. y, jplOT 

t SN 

0C3F 

2 F-jrWAT (215) 


jnN 0039 
I f N 0^4 1 
If N 004? 
X SN 

ISN 00* 3 
I £n C 04 
! SN 0C«*7 
I 3N O C«»f 
I fN 0:49 
I Sn C C 5 C 
r^N— 

!.:N c: = i 
I UN OOjl 
liN 0 0*: 4 
T “5N 0C5* 
I UN OC*>c 
ISN 0C57 


JTMO'>V=«, CAiL THCUfiV, PFAC IM CATA POINTS, OC NECOe«''iCN ANALYSIS 
JTMONY^), C*LL ■'HtC'PV t'NLY 
5r 0*’’‘ ''0 WC-.BSS-.ION ANALYSIS AtThCUT CALLING TNECOY 

1/ JPL_T»l, TM*"N A iAwPLF CUtVf OY Tpr jfp «jth SOLID t ABTM LUYE NUNBESiKK A 
/c»^0 tAT, ANGLi Ib ^^CTi'f-C rtV THFHhV 


IS ^^CTi'f-C rtV THFn>^V 
•£C. 2) GO TO 6 
• cO« I) iC Tc :» 


IT ( jTMiipr 
4r ( J’«’HOWY 
GO TO 4 
CALL 1 Ht J-<Y( JFL JT t MT ) 

CALC LAiJTA 

CALL Cn»»we«( Tr-TAP T.t cr^n.OT.MFXP I 

CALI 
njU to t 

call ’^HFUKVI JFLOTi MT ) 

GJ T J f 
rONTINCT 
CALL LArjTA 

CALL CCrP£S( TSTART •T£HL' • OT t MEXP I 
CALc Wf.G«>£?(MT| 

CUNT XNLL 
*>Tl)P 


I SN 0 00 2 


! 0 00 3 

ISN 0004 


I *N 000 * 

ISN o:-cc 


ISN 000 7 
ISN poce 
I SN 0CO9 

'rsTr-coro 

1 SN 00X1 
ISN 001 c 
ISN 0013 


ItN 0014 
I.SN 0 015 


X5N OCiC 
ISN C 0 1 7 

? SN oci e 

I iN 001 o 


.vUMwCUriNb AQC4 1 T > 

THn:-^ri>R3ir-TNc cnMncTE^ rwr CC.TTNC op the arg«iment cp tip <i.e. it 

COMPJTL.i THL COSIN'" Or (L - «( APGUHLNt OF PER IGE Sf 4 e Ar, ASCMAiyi 4 C» (mean 

anomaly) ♦ M'CNJDL QF 01Sru»BlNG BODY) - M-* ( NODE CF SATiLLlTLJ ♦ TICAL LAC 
AN OL C ) • 

IMPL ICXl R^AL••HA-H,C-^) 

C IM 4CN/HL<A/CP JLNM.cPSLNJ* ACGII4 , AR GU S • AP GOT 4 * APC f 
aiNMJN/tiLND/«LAG.'iiOMc.GAM.O£LHN#CM£GAS*OeLSUN.5.5*X!0*A,£.Kl.*NOOE. 

1 IZ^PCj**! >^01' N»XI 5 iN*TMS#irNOM*«NOn#F# Y’4E ANM tX4EAN‘^*K40''7v,KMDCT<« 

2 tiMGOCT.McXO 
OU 1 C LM«l t2 
VL-^sLM 

00 1 0 LPsl , i 

VLP*LP“- I 

OLj J 0 L ’O*. • ? 

VLQ*LG - 2 

A?*=(2.0l0 - 2 ♦ :00* YLP) 4 < XNDM*< TT-Tm 3 )40£LMN) ♦ VLQ • ( X VC OT M4 ( TT- 

1 TM j ) ♦XMt.ANM I ♦ VL4« ( CMGOtjT* ( TT-TMRI ♦GNFGAM) - YLM«( Q»( TT>TZ£P0) 4 

2 XNODc) * LPSLHM<LM,lP *£0) 

APGUM( lm,lM*L0) =DCOS( A2) 

A3s( 2»0r 0 - 2 • OuO* YLP) •( XN05«< TT-TM5 14-DELri’M 
T*T3 )4Wu AN j ) * VLM -omegas - VL M* ( rj • ( Tf-T 2F PQ ) 

2 LPiLNSlLM,LP . lOI 
An GU 5 ( LM • L P • L C I s OC CS ( A 3 ) 

0 Cur;TlNUS 
RETURN 
tNO 


VLO ' ( » vocji c*< T7« 
*NC0r) ♦ 


X SN C 0O2 


ISN 0C03 
TUN 0004 


ISN 0CC7 

1 SN 0 00 8 

2 SN 0009 
ISN OCIO 
I SN 001 1 

-TSN “CTJT2“ 
ISN 0 012 

! SN OCI 4 
ISN 0015 
ISN OOl« 
ISN CCl 7 


SOB«-OUTIN£ AfiGDOT 

THIS 3UnRnUTIN£ COMPUTES THE AKPrliLAP SPEED OF EACH CCNFTITUENT II. F. IS 
Tut TiMc DERIVATIVE OF THE ARGUMENT COMOUTEO IN SOBRC'JTINL ARr,>, 

implicit «lAL»'>( A-P,C-^) 

CQmM JN /ULKA/EF SlNM •£ P5LNS • AR GU4 # ARGL S • AP 00*^4 • ARC CTS 

<L*f-*Q.nMF-.AM,’3ELMK,CMECAS,OELSUN,R£*XT0. A.t ,«I,XNOOE, 

* TZCnC • X I MOON •XISUN*TM3«XNDM«XKDS*F*XMEAKM.X4EANS*Y4CC7 v, >f MOCTS# 

2 UMODOTtMEXP uuia* 

DO 1 C LM=1 ..2 
VL 4-sLM 

00 1 C LP»l • 3 
YLP=LP - 1 

DC 1 C LO*l , 3 
YLQsLD - 2 

APGD TM(LM* lP *L 0 ) » I 2 . UDC - 2 • 00 C * YL P ) • X NO M ♦ YL0*XMCCTM - YLM«Q 
1 ♦ YLM»OMGDaT 

0 Cunt INlJt "^* *~^ **" ^ * * * ^ " ^^^ ~ 2*0P0*YLP)*XNDS ♦ YLC*XMCCTS - VLM*C 
RETURN 
END 
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I5K 




iU.«**rM'T2r<r AC'5L?r<TT> 


c TMir il;’' ► I ^ Thc I NL ^^0 CIjInE TF (L - T • GUNFM Vf 

C AN %NnMALV| ♦ O^CMTAN AlirMAl.V) ♦ M*(f4CCC Cf fTS^UrPlN'j HCDV) - 

c •-•CNnut. CF pjr u>r *n “LjcjuTira cata. 

c 


c 


I rM 

Z SN 
I ON 

0002 

COCA 

ooos 

IMPL KI^ H.AL**(A-h,C-2» 

./IMF Nb HiN C>W(&f2« }>f-.NM(2*3*2l*C5'5(2#3*2>»SNS<2*3*.') 

Lf'N tCN/lJLM / XL AO* U AM* r^LMN* CMC GAS# -jElIUN* i;#*Xl':#A,> •xl* XNCOfc 

1 TZl PI *^ i MO *N. X 2 S INf'^M i, XNPM *XN01- ,F * MM^ANM, XM£ ANS . XAC*"TM , wM^rTS, 


* 

2 aMGorr.MEXP 

: *N 

uoot 

COM «r»-/MLXr/Cl'M* *NM, ct at 3N." 

I ,N 

ccr'T 

OO 12 LM*1 *J 

I Sn 

0 : 0 £ 

VL 4=LM 

! SN 

0 OO'i 

3U * C LP«1 • r 

1 SN 

OCl 0 

YLP*LH - 1 

I £N 

001 1 

JL 12 LO* 4*2 

UN 

OOi 2 

VLUaLC - 2 

T-rv- 

C0l~!' 

A^^C2.■3^•■3 - 2.0C?* VLP|T( YN-JMrfTT-TUGUDELMM A VLO • ( v vlC’ m# ( TT- 
* TMG >♦ anm 1 ♦VLM* ( CMOC-JT'- f TT-Tmc) ♦ lMt'“'AM» - YLM- ( 0-< TT-rZcRO) ♦ 

2 XNJDw) 

I SN 

uOl A 

CiMCLM.LP.LO )=r>COS< A2> 

1 FN 

C '7X 5 

SNM(lM*LP*LO )*:iblN( Ai» 

1 LN 

0016 

AS>(c*OuO - 2.0n0«vi.o»*(XNr)0*<T’'-TM'C UJELSHN) ♦ VL 0 • ( Xm COT 5v < r T- 
1 TPS UX./EAN i »♦ YLM^OMFGAS - YL M h < 'a • ( TT- T 2E PO ) ♦ XNCCtl 

7. 

OCl 7 

Cbi(LM*LP*LO)=>CCS(AU 

TTT7- 

7TTTF” 

5Hr(LM*LP*LO )=OETNC A2) 

I £N 

: n c 

10 CON^'INUL 

t$N 

oc »0 

PC TURN 

I £N 

0 32 1 

ENP 


t SN 

00/£ 

C 

SUrtP J*>TI Nc. tJLM 


■ - 

c 

c 

c 

THTT •UORaU^lNI COMPUTES THC :*C^oniAL PART CP The TIP, 

1 SN 

000 1 

c 

implicit PCAL»9 (A-H,n-ZI 

I SN 

OOOA 


OIMLNS U-N Ei< 21 • XMCOM ^* X, 31 ,x SUN ( 2* 2 * 31 *Y*CM< 2 *3* i ) . v»c b ( 2* 3 • 2 

I ?N 

0005 


Ci*MM 1F*/BL<G/H: *XK,CA'NTM*CFNTS«*?HM* XMCON* XSUN*YXW ,YX ;.* XKf» 

I ON 

ooc e 


HI ( 1 i-i • :/: .ooc 

; SN 

0C07 


iU(2Ul. 30C/U.0no 

rSTT 

“OrTJH- 


— PETUmJ 

1 SN 

C03S 


cnO 


I <^N 0002 


I£N OOOJ 
ISN OOOo 
I£N OCOS 
I £:n OOOfc 


liiN 000 7 
I«N COCA 
ISN OOC9 
I SN OOi c 
t SN OOl 1 


lUP Ot»T i Nf C Cr F F 

c 

c This £imt.'£j*jT I nc cnMotTes thf part hf thf tip i*hich maj ^he for** 

C Cc N5TANT/( A*^* 7/3 ) * 

c 

c 

IMPLICIT REAL*^ (A-H*0«7) 

UiMENt It.N 91 ( ii > • KMOONC 2 t J • 1) . ''S'JNf 2, 3 • 3» #VKH< 2t.3«3)fY»(?(3»3»3) 
CUMMOFi/iiLNH/eit AK.CFNTPfCFNTSfS'JP* XMOHN • Xf-H A • YK M , V K £ ■ * K P 
COM-AON/eLKO/ XL AG*Q •TMtGAM, OELMN* CMEGA5>»OELSUNtSS ♦Xl7tA,6,XI,yNr.3Ef 

1 T7£Rr . xihCOn. <! Sun*Tms* xnom .xno «. f* x^t anm, xpE AS bt x "CC “v. ■ mdcTs* 

2 OMGD*<T#MfcXP 

A^sOSCPiMAJ • 

CFNTMsrtJ 3.22 3200/( <AA*31«A2) 

LFNTb«i72.571^D0/( ( A#*3I*A2) 

HfcTURN 

cNt> 


ISN 0C02 


I3N 0002 
! FN c oo A 
! £N 3 0(rr 
TEN OOCF 
UN 0 00 7 
Ur4 OCC^ 
T«N OOOG 
1 SN 0011 
UN 0CI2 
UN 0 013 
rSN- 0 C T T 
ISN 00*7 
I 3N OCl e 
ISN CCIP 
ISN 0 0£C 


SIIBFCUTINE CCPPESC T ST APT • TPNO • OT , M£XP| 

C 

c THr 3 ‘ ruripn'jTitir- IS U5£C cnlv nv *5U0 RCJtine plqtcr. p’jTj in the values 

C OF APnAY F'TPACK(1I, WHICH OETEPMIFJES WHETHFC OP NOT AN E» t>£ P I m tNT AL DATA OO TM 
C SHOULD 9E PPINTtC AT TIME T< I J ( NTP ACK ( N ) =TME NUMUEF CF TmE LXPERlHENTAL 

C CATA POINT iF Yt3« 7ERC IF N»J« ) 

C 

c 


implicit PEAL*A<A-h,C-ZI 

01 ME n: ION "^EXPI^OO ) • XI £XP( 200 I • T(70C UXI NCL (700»*NTiiACK(700» 
•OT8WnK7BLKC/Tr<P"fXTESP *T* X^INCL .NTPaCK 
MTs( TtNf -rSTAt-T )/0T ♦ UIDO 


0 J 10 I = I • M-T 
X 0 N7L ACK < I >=0 

IF (Mtxr .EQ. *»> GO TO 12 
UO 11 l*l.MEXP 

Ns< texP( I )-T ST API I /r>7 ♦ l.SOO 

IF < N cLT • 1 ) N = 1 

IF IN ‘.GT, MT1 N*MT 

11 N*iPACK(Nlsl 

12 CONTIfjUL 
PtTUPM 
tNC 


47 



»5N OOCe 


I t.N w ,* 
I sN 0004 
! 000? 

in* 

!i.N 0 00 7 
I iN 0 : c c 

? L.N 0^0*i 
I SN c : 1 0 
1 CN 00 1 } 
! -:N 0 01 * 
! j>*W 0 010 


Mrt*'t3UT INc. ) 4 <mm ,01 

THi^ zcnr^nijT iNt* is u«lo to ^r^'O ts '^ata fo»- lse in :»*l v.il^iplf 

ANAIVS;'- CAhT P^o OUT” * •.tJ-IDU^ ! NE PECMES* 

AL C A-M, C- i » 

•^rAC»A 

Ol Mt N: ion t.PSLN'* I ».' • t# 3 I # EFT'E N** ( 0 • 3 • I • APGli*'^ «i t 3 • J ) • ACf,t ■, ( ^, ), ) ) , 

1 A».GC1M( J* 1* 3 ) tA»iGJTS< 2. ?) 

nxMrK:To>< pi (2 > • irvcrN(* .j# n • xr*jN( 2t3 # j )#▼*«•»< 2 * j.: > ttr' ( 2 . i*o» 

01 MlN.- I( 4 T*- »I*»C .0 ■' » # « :i*P< ?0 0) • 7( 700 » » »I NCL < 70 0 » .N^K ACr ( 7001 
L'IMlM ION C( J. 2) .Ai( 2* 2) »e(2# J I .CLPOMI 0,2 3 #GLPQS (3*^3 
LI wTAi :r.r< 2. 3* 3 3 •'jNM( 2 # 7 ,0 3 #r$G( 2* 3*3 3 .SKSC » • 3* i 3 

Cl'PMUN/Ul-X A/PP«^LN'^ 9f P<-LNS * AOGl'M , AOGUS • A»»G07 *^ * AfiCOT? 

CL»MMON/eL»^ Rl • AK *CFNTM,CPNT«,G*J3^* XMOON* Al tIN* VrM.YIt . 

COM^ON/MLKC/ AP , Al M»‘* *T ,X I NCL* S7f3ACK 

CU^’M»)N/('Li<.f«/ ■LAf.*O*C4EGAM*0tL»^N*CAEjAS* OEL #S5**XlO*A*f, ,Rl»XNCDEt 

2 T7rpr ,X!MCON* XI rUN.TMC* XN'^P#* no *“• F* XMEANM. xml an: • X von v, wnr^CTS, 

J C wSL:»T *^<» XP 


1 tN 

00 1 A 


C j( ►* JN/lti.K.r/Cf V .INV, CS3, 5Ni 

1 SN 

wOl ? 


CUP»^CN/NLN 0/ c, a: , JL’QV.CLPG*”- • I NDF X 



C INO 

“X Iw TL. X=FP 74ACN OF 7XmE 

t SN 

001^ 


IN'JF ■* X ♦ 1 

1 SN 

001 ^ 


TTsT£»rp( iNDt X 3 

1 SN 

00 1 c 


CALL APGl.5T(T“3 



C D1 

Tr»E; LUNAf •IDES. "F LU«® ThF LUM^GLAR TICES i 

I SN 

001 c 


1-0 

:SN 

co.:c 


DO 73 lm = X 

I SN 

coil 


YLMslM 

I SN 

0021 


OL "70 LPsl.r 

1 SN 

002 > 


f)t 70 LO* 1 » i 

1 £N 

002 4 


U I « 1 

I SN 

Ovj-« 


f>( I 3 aCPNTM«Hl fLV)*C(LM,i).Al(LV,LP3*GLPOV(Ln 




1 /■{ AK SC7M(LP# LF *L 0) X ? ^ J 

I sN 

0C2t 


l = I ♦ ! 

1 SN 

C027 


J( I 3=CFNT^*L)X(L 0*C(UM,23^AI ( L P . Lf* 3 « SLPOM ( LP 




1 / ( Al« Gl TM( L V ,LP *LO) T ® 3 

1 IN 

ocie 

70 

CUN"^ |MJ( 



C tlO 

The HU. AP TIOEt 

1 SN 

c 


Oo 7 I LVsl ,L 

I SN 

CC3C 


YL -^sLV 

Tsr 

0(771 


OD Tl LPsXf 7 

1 ^N 

30 U 


CO 71 LO=l .7 

I SN 

00 3 3 


IF (LP .uN. 2 .ANU. LQ .EO. 73 GC TO ?X 

: DN 

CO jf 


is t ♦ 1 

: SN 

003t 


L»( I 3*C ENTS »U 1 ( LM ) • C( LM, 2 3 •P'CLM.lP) •'■-LPOt (LP,i 




1 ( Ai5 GO r ■ ( 4-V, LP.L03 *1-? 3 

1 5.N 

0C3 7 


is I ♦ 1 

: SN 

CC 3S 


I 3 =C FmTI 'C 1 (L-4 3 C (L«,2 3 • H ( L« • LP 3 vOL^ 05 ( LP • i 

— 



1 ( AF or TTt l*7,lp*i nx ■rsj 

I SN 

0 JIG 

71 

ct;M iNUc 



C PM 

IN C21LRV10 INCLINATION 

I SN 

o:4o 


1= 1 ♦ 1 

I SN 

c :a 1 


u( I ) = »l-.XP(lNi;Ex 3 

1 CN 

004c 


Rfe TURN 

I 3N 

0 04 7 


cND 


I SN 0 002 


X 3N 0 00 3 
OOOv 
I SN 000? 
I SN 0 C 0 c 
TCNonj 

f«N C0C3 
ISN Owli-' 
1 SN O C 1 1 
ISN 0012 
ISN 0013 
ISN C 01 4 
ISN 0012 
g 


I TN on* 

ISN o-oie 

1 SN 001 G 
!SN 002C 
v»02t 
ISN 0 0? 
ItN 0 02 2 
TSW“T72* 
ISN 0 02 = 
ISN 0C2e 
ISN CC27 


SUni^04»TlMfc cCC^UN 

THT“ rVDPr»UTlNF CnWPC^ES THT rCCENTPICITr FU^CTlCNS F2F Thc PCCN AND SUN. 
Implicit realms < a-m. c- 1 3 

DIm:mICN C(«*3>*AI(2*33*E(2.7) *GLP0M( 3,3) •GLFOS<3 • 3 3 

COMM.jN/ftLKG/C . Al ,H*r,LOQM,GLPOE*INOEX 

£MCON = O.02A‘iC0 

^A.PQV( X • 1 3 s( - C .*5100 3 «E^ODN 

GLPQMl 1 , l!)=l .030 

GLPOMC 1, * I =( 2 . 50 0 3 »SMDON 

CiLf' JMl a « 1 ) 3( !• OCOI *£4CCN 

GLPO >(( «;• 23 «l . C“»0 

uLP<JM( 3) s( I .SCO) •FVCCN 

C«.POM< 3, X 3 *( 2 .SOO 3 *tMUON 

OLPOM( 3* 2 3 *1 .00 

OL PUM( 3, J) s( -0.500 3 • EMCON 

.c 167D0 

OuP'3 j( 1 • 1 ) *( -O.SDO 3 *eSUN 
GLPUS( 1 • 2) «X .030 
GLP03( I . !3 = ( 1.30 0) »ESUN 
uLPUi(2.X3=( l.GCC) *t5UN 
Gl P*l3( 2* 4 )»I .COO 
ULO)S(2* J3=( X.SOO)*cSUN 
GlP03( 3* 1 ) = ( 2 .SCO 3*' ESUN 
GlP X5< 3.2>=1 .OOC 
glpos( 3 3 = < -0. SCO >*esuN 

kETURN 

tNO 
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1 


1 


I ^ 0002 


30HttGUTlN£ !NCL<«l»F) 

THUTTTJPtlOUTlNr CGMPUTcS THe‘<:5CCN0 DEGREE 


: FN 

C003 

implicit r 

I SN 

0004 

OIMCNSlUN r 

ISN 

0<30* 

S*Oa INI XI > 

ISN 

0 J06 

C»OCO‘ ( XI 1 

1 SN 

0CC7 

S2=D*»2 

I SN” 

C mC H 

C2=C'* T 

I «N 

0039 

F ( 1 • 1 1 = ( 3 • 



F( 1 *21 = (-3 

I SN 

0 311 

FI 1 . 3)=« -3 

irs 

001 z 

F(2.n-{ 3. 

1 SN 

OOX 3 

FI 2* i»«I 3, 

1 SN 

00X4 

FI 2* 3)s( J* 

I SN 

on 5 

HfcTUPN 

i ^ rv, COIL 

END 


INCLINA** !TN F’JNCTIONf. 


1 iN 000? 


I £N C 3 ; .■* 
I fN 000** 


~t; n ' J TO? 
I SN OCCt 
X£n 0007 
IbN oooe 
I SN 000 9 
!*N 001 C 
X SN COl 1 


3UHK(X»T|Nt LAO I 

TrtI S EU0P3UT !KE WcAC3 !N ThF 
rnjA DLGKtEb ’'•J BA"«IAN3« 


SCLlD earth lag angle XLA'*. ANC CONVERTS IT 


CuMMON/ML<D/ XL AG . 0 tO we G A M , OE LMN • C M2G AG * OEL 5 UN .SS.XICtA.E.Xl, *NCOE • 
1 T 7.i.HC • X IMOJN.XI SUN.TM5* XNDV »X NO 5 • # XM^ X N M ♦ X »#E AN :> t X MCCT s* * X M JC T c , 

i UMGOQT, MEXP 
RF AT* ~T3 • 1 ) XLAG 
FUPM AT < <3K 1 0 •£ ) 

FJHiAT^?////////.23X,*50LlP EAH'-H LAG ANGL t = • • F fi . 3 • 1 X , • DEG9E K S • ) 

XLAG=XLAO* F 

OtTtjRN 

cNO 


I SN 000 2 


TEN 0033 
I ?N 000- 


:*N 


I SN OC37 
! SN : 0 3 ^ 
!Fn 0039 
I£N 0 310 
I «N 0 011 
—fsK TOT'S' 
I«N 0013 
I 3N 001 X 


FUHhOUTInL LAG2 

THT5 TUSnrVlTlNr 3ET3 THE Ft^roUENCV DEPENDENT TIDAL LAG ANGLES INITIALLY 
.DUAL TD 7CF-3. 

IMPLiCIl RFAL*®(A — H«C— Z) 

OlMdN*- ICN £P«LNM ( 2 »3, 31 • EP9L ( 1 * ^ • APjUMf 2, 3 »D> , A^GUSI 3* 3« 3) • 

I APGCTMI 3* 3» 3 ) » AHOOT S( 2» 3* X ) 

C )MMCN/HLK A/CFiLNM .FFSLNS. APG’JMf ARGUS* AM COTP . A RGC T ? 

CC>‘f'raT^/l*LisD/ XL AG,Q ,n M“G AM , f^EL^N* TMCG Ai • 0E._ SUN # :*F *X J . * • E* * 1 • X^CDE * 

1 r/LPL* X IMJUN, Xi bUN*TMC- • XNDM*XND “*F, XM£ A N.A , XME ANi , X V LC . M, XMOCT.c* 

2 tJMGDUT.MEXP 
OJ 10 LM*s 1 , ** 

OU i 0 LPs I • ^ 

Du \ 0 LU=l *3 
EP?LNM(L ULP *L0) =0 *000 
FP 3LNi (LM, lP *L0) =0 . 3CC 

10 CDNTINUw 
MuT'JMN 
END 


OEIGTNMIP^^ 

OF POOR SOAI"” 
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I 


loN 000 ^ 


75 t^a' 5 » 3 r 
I :: N Q OC •* 
lit* 0 JO* 
I oN OOOt 


I SN 0 “'ey 
I£N 0008 

I SN 0 OJ C 

: CN oo;i 

I SN 0 31 c 
ISn 0 313 
IsN 001*» 

» 0 31-3 

1 ?N 301 f- 
T' 3 ^ 07717 ' " 
t ?N OCl e 
liiN 00«9 
! SN C0?C 

I iN ooai 

! SN 005 £ 

I SN 0 3.:? 

! SN OOJ A 
ISN 0'*>?5 
I SN 002f 
IcN CO-7 
1 3N 0 02 S 
ISN C02-5 
IfN 0o.?0 


I SN 0 0 3 2 


ISN 0303 
I SN 0 o r 
1 SN :oos 
I 3N OOOc 
T SN OCO ' 
ISN onot 

UN co:^ 

I Sr< 001 C 
I SN O0» 1 


5UHf<CUTlNc L^SOTA 
C 

c 

C THmr01'.l>^OTlNp~PerA0S Xn, ThCN PKIMS out. T»-E □e«*ERVCD tip. XlEJrp(j) is 
C THc. fJdSEl*VEO TIP (»N AhCSECONOS) AT T T ME TEXP(J) (IN DAYS). 

C TOAY IS A COURtCnON FACTOR* WHICH MAKES TEXFCJI START AT TFE tiECIKNING OF 

•C The YtAR {JAN 1 AT 0 HRS UT = 1 ). TCAY IS NECFSSAPY SiNCf OFTENTIMES TPE 
C TiMc fir EACH L'ATA POINT 15. GIVEN IN OAV3 AF TEP THE FIP£T CnSEPVATICN (l*F. 

C TEKP(ll«0.0 ). 

c 

c 

IMPLR'IT REAL* 3(A-H, C-2 ) 

01 MEN. lUN TI XP(200 > . XIEXP(20 3I . T(70C ) . KI NCL I 700 ) *NTR ACM 700) 

CU.M.MON/CI.KC/ TEXP *XlbXI> *T#Xl NCL* NT RACK 

CUMMON/tiLKO/ >«LAr,*0,UMtGAM,CELMN,CMECAS*OeL?UN*SS*XlC.A,F.Xl,XNCOF* 

I TZcPl. . » I MOON • XI Sun. TMS . XNOM ,xNOS. F. XMEANM ,XM£ ANS. X moot m, X moots* 
i UMGOOT.MEXP 
PEAO (3.31) mfxp.TCAY 
3J FQOMAT ( 16 .F 1C. 3 ) 

" ■ DO 32 J = 1 . ME XP 

REAO {‘S.3'4) XIEXP( j) ,TExP( J) 

TLXJ* ( J IsrEXri { J ) ♦ TCAY 

12 CONTlM.rc. 

34 format ( i0l4 .F ) 

WRITE 

^ FUPMAT ( IHI ) 

*H ITE ( f; . 2 ) 

3 ■■ rUPMAT ( ///// .“OX . •5ATELL I 'E TRACKING CATA*.////| 

WRITE (e.lE) MEXP.TOAY 

13 Fjj.RMAT (4 1 X, • ( • • li> , • DATA OCINTS. TO A Y= • , Fa . A , • CAYS)**////) 
wPITL (f>.A| 

rOPMAT { 4rtX, *T I N*L • , 1 2X. • INCL INAT XON‘ ,/ ) 
wo I*^c {-5.3) 

D FORM AT C *7x. • ( C AYS ) • , I 4X* • ( AOCSCC) • ) 

WP ITL ( O. * f ) 

I^ format (3aX,*(JAN I, 0 MRS UT = I)*.//) 

OJ 10 J» 1 . Mr XP 

10 WRITE (fc.l*) J*TEXP( J) .XIEXPI J) 
i 1 FIJPM AT { 33X.If.SX.FlC.4.l0X,F10.Aj 
Pc TU«N 
ENO 


SUHRfM*TfN- LOVr 


C TmI3 IUl'POLTInF RcAOS Tn Tnr SECCNO OFT, REE LCVE ANC LOAD 

C Ot.FORMATION COEFFICIENT CF THE ''iflLlD EARTH. 


implicit RLAL^oIA-h.C-Z) 

UIMENS tON bl (2 » . XMQONC 2.1.3) .x£UN(2.3 .3) .YKM(2.3.3 ).YkS( 2.3.3) 
C')^*MO^;/bLK n/bl . XK .CFNTM.CFN* 3. SUM, XMCON.xSUh ,YKM ,YX S, >KP 
PcAU (S . I ) XX ,XKF 
- I FORMAT (3r 10.») 

wr ITd < 0.2 ) XK , XKP 

2 FUn.MAT (//////// .23X,«S0L10 EARTH LOVE NUMHER^* • F a . ** , l OX , • LC AO OEF 
lURMATiCN CTEFF I CIENT* • ,Fa. 1 ) 

•return 

END 
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1 


i 


iUH«n‘/TlNh LLVMUMIMTI 


T SN TTCT3 2 
1 SN 000*4 
15N 0CJ5 


C THI.1 CLTmTjTlNr T4<F5 TmF COiTrF t <" 1 HnT J ANO ▼h£m FCLNC IN 

C sUnWijUT i:.C and -*CLVE3 for ^Ht FR*JO’jrNC>< Oef»eNPEKT LOVF NUMHEPfi 

C UAo ANvjLfrSi ANU Thti ^ CRWORS. !T ALSO FINDS THE HCE AN TICE P A R A#«£ T E R£ , 
C AaoUMINi THf. JULIJ CARTh LCjVh NNMBET AND LAG ANCLE ALL *fSC<»N. 

C««*#*Al.L Tl.'fe.S ARt AS5UMF0 Tj yt ONC rEG^'Et. ONL^* HtNCu ALL CThfclP TIDES 
C AdbilRtiEO IN^n T»-E <iNO OEGHEE ■*'IOE«. 

C 

c 


1‘APLiri- r^cAL* “( a-h,o- ») 

PcA;.**. VV( *^0 ) . Vc ( 70 I * ANSI 1 0) 

ji mln: xln lpsl -imi 4i, j . ji .e p 3lnc < . 
i AWGCTM( », 7, 1 ) , ALC.OTSI a, 5 • 7 ) 


• J • 3 ) • ARGIJMC 2 ,3 • J)«APoLX«(Jt3«3)f 


I ?N 

OOOt 


CIMr NV 1 LN E* 1 ( » ) , XMCCN ( ? • 1* 3) . XSUM 2.1t3» .VKM(a. I.Jl*V»<S(a*3.3) 

I SN 

GOO 7 


.01 MENS I ( H Tl 20 J ) t XI f.XPIcCOl t T<7 OT ) , «I NCL (7 30 ) • NT PACK (770) . 

I it< 

OCOP 


J iMuLL U.N C( ' t 7 J . A X < 3* 3) # H(2 *3 ) . ULPQM ( 7 • 7 ) # GLPQS ( 3* 7 ) 

I £N 

000^ 


DiM-N^'lLN NN(70)il3AVt.(^0> 

TEN 

crro ■ ■■ 


DIMrKriDN T!0E(2.3*3) 

! jN 

('0*1 


DiMEN.'. K N SAT( 14 » 

J 3N 

0012 


COM ^C^J/HLN A/tP3L 4 A P 5L NS • AP GUM . AR G'JS • APGOT M • AR 0 C T S 

I iN 

001 •» 


C’JF‘^Or.’/i<LN|j/ lii • XK, C^NTM ,CFN“ “t , S'lM , XMOQN, X5UN • VKM , YK « , > HP 

T SN 

OCl A 


Cl WMON/ULNC/ TE XP » X I F XP , T , * I NCL # NTR ACK 

T Cf«4 

001 5 


CUMyON/HLKO/ XL AG ,0 .C G AM , OcLMN • C M 3G AS • P cL bt »N , SS tXlO,AtF#xI,xNCOE» 




1 V2E R: V* X I MP.-.N, X 1 =UN* TMS*XNDM,XNDS,F,XMrANM,XMcANSfXMr^rTM» XMOOT*" t 




2 OMGCi’T* Mf. XP 

TEN 

UDTE 

" 

- CrJMMDN/nL* j/c. A1 ,U .OLPO*# •GLP0‘‘.* INDEX 

I SN 

OOi 7 


C( PMCN/ULN*-/ W • VL i ANS* I b AVE, M, K 

I «N 

001 ? 


cjHMCN/uLN :/ 3A r 

I SN 

C •- 1 P 


CATA PLANK/4H / 

I £N 

0 02 0 


JATA “Mi 1 2 , TMi 72, r M2 12 • TM2?2/4hri ,4HKl ,4PM2 ,*HKi / 

I SN 

00c 1 


DATA YSl 1 2 *T SI ?2 *T453 1? fT-aa ?i- /4Nr»l , 4 mx 1 S ,4HS? / 

I SN 

0 Id'k 


PI = i. 1 41 bP2653“15D0 

I SN 

002 3 


TL A ‘. = XLAG/r 



C 3 

rr THE ACDITIVC CCnSTAnT in Th-: T!P “DUAL TO xio 

I SN 

C 02A 


X£0=ANC( 1 ) 

I £N 

ooas 


•FITE {<- .ibl 

! 5n 

core 

1 i 

rUPMAl ( IH 1 ) 



C p^*nF *;u1 TMu headings 

I SN 

0027 


•RITE <3iH1) 

! SN 

J02 d 

o 1 

i ••-tA^-uCcAN TXOt PAPAMt TERS ) 

tyTT 

'OOSO 


*’RITF IF# 32) 

ISN 

OOTO 

H2 

FuPmaT ( iO X, • I as SUMI no 2ND DEGOCr ’’IOCS CNL v » • 1 7 x , • ( A^DI 'M ! NG 2ND 
TiJcCM 

1 SN 

0C3 1 


•RiTE (6.13) 

ISN 

00 Jc 

H J 

fopmat ( 7c#h, •lnly and solid earth LJVr*) 

ISN 

0033 


• LITc xk 

tSN 

OOIA 

34 

FuPmaT ( ax , •NJMUcR k2= • *F5* j * 1 X , ' and OOLID*) 

T SN C 0 J5 
'Ttw- • • 

53 

• RITE (c,,3E) VLAG 

FORMAT r-x-jx, XEARTH TIDAL L AG= ’ t F-I . 1 , 1 X , • DFG ) • • / ) 

T SN 

0 0 3 7 


*R *TL ( 4^ • 40) 

ISN 

00 36 

•*0 

format (IX,*L •« P ij TIDE 02G*b0DV LJVE NUMBER STL. c-HKCR LAG 
IAnGLL 5TD* fPROP C ‘TD. COflOR PHASE *=TC. Epppp 

I SN 

0019 


*Tc (f>,39) 

! SN 

0 o**c 

39 

H JMM AT 1 31 X, • (DEGREES) • , IX • • ( DE»iFEEl )**7X,»(CM)«,4X,MCM)4 ,ik, 
1 • ( : 6GPLt_S ) • • 4X , • ( DEGREES > • . / ) 



CT 3ET NN(JJ so fop tidcS not SCLVED iiv THE MULTIPLE LtGPtS«ICN IN 



C o 

■JtJR.JUTiNc HCOf-Ec* AND NN<J)al FPP ThCOE VMICM APE SOLVFL' FCR 

ISN 
T SN 

0C*4 1 



0 0A2 


DO 3 0 J = 1 • MV 

I SN 

00A3 

2 0 

NS( J ) = 0 

I SN 

COAh 


OO 21 J*I *K 

I SN 

OOiXS 


IoS= H AVE ( J ) 

» SN 

004E 

El 

NN( I S*)=l 



c # »#: *#:LI)NAR an.) LU N I “*UL A P TIDES* »<•' 

c 

C PUT NAME OF tidal const I*'’U£NT IN APPAV TIDE 

I SN 

0047 


00 d9 LMsl ,2 

I RN 

oo4e 


DO AV LP* 1 • 1 

I SN 

0':49 


00 39 1 0 = 1 ,3 

I SN 

1050 


7IPE (LM*LP,L C» =MLANX 

i SP4 

"ITOSl 


TIDt C . 1 • 2) =TMI 1 2 

1 SN 

OCSc 


TIDc ( 1 * 2 . 2 )=TM 122 

XSN 

0052 


TIOEIc.l, 2)=TM212 


P^G^i ^ 
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TICAL 

ARE 



: ?.N 


! aN 

OCaE 

I SM 

00^*3 

I SN 

::s7 

1 SN 

DCad 

I rw 

■J . a S 

I SN 

DC 5C 

ISN 

DC 3l 

I ‘.N 

00^2 

I LN 

OOe 2 

I ?N 

OGha 

1 SN 

OC -SE 

I SN 

3 0f:t 

I "N '006*7 

I £N 

CO'' c 

I SN 

0Cc9 

I SN 

0 : ■»c 

I SN 

j:7i 

1 GN 

DC 7 ? 

I SN 

o;7L 

I SN 

DC77 

I £N 

007E 

1 KN 

0074 


ZZ^O 

I SN 

0091 

I LN 

o:a2 

I SN 

oo*-* ’ 

! SN 

Coda 

I SN 

OCAS 

1 SN 

OC AC 

! aN 

3 •: *> 7 

TSH 

0383 

I SN 

00H9 

I SN 

3C9C 

T CM 

0091 

I £N 

C 002 

ISN 

009 3 

I SN 

0C94 

I YN 
ISN 

CC95” 

0C9o 

ISN 

0 :97 

I SN 

ccie 

! ?N 

0099 


! SK C I 


0 


I£N 0101 
1 »N O I O 
ISN C 10 J 


I 3N Cl CS 
ISN 0107 
ISN OlOd 
ISN 0 1 1 0 
X SN out 

ISN 011^ 
I SN ~jn 7 
ISN 0114 


V I DL (*•_•-) = 
LLL =2 
1 = 0 

OJ 6 0 L'^^i t 2 

on 'OC LP=l .3 

VL^^SL*-* - i 
- ; 

L\J 6 0 LOS 1 , 3 
LI.OsLl - 2 


I»I ♦ 1 
J1 =VV< I 1 
clsvt( I ) 
risT 

1= I < I 
Uis-VV( I ) 

I <£ =vr ( J ) 

C iH JNLV *' Ht. i.rf£C’’IVL LUVL NU^flCP !T SOLVfTD F CP * P»Jt THE t.»-FECTXyF LAC ANGLF 
C ANJ tP"iJP tCUAL TO ZEPU 

Ir fNN(i) •CiC. 0) U2=0*CD0 
IF 1 NM 1 ) 0 1 Eis?.COO 

c IF THIS -»A? nct scuveo nt , go tc 02 

IF INf lM) •&€. 0 .AND. NN( I ) .CO* 0) 0(J TC 62 

UJ = 'J 1 *‘^ ♦ U 2*»2 

C FiN;^ THE i.fFi.C“lv£ FHEOU6NCV-or:PFN'Ofc.NT LOVE NUMe£P 
YKV(L.y.LP.LO)=OSOPT(g?) 

-J.*=U 2 /L 1 

C FlNt> THE EFF'JC'^lVE F ^<cOUr. NC Y-OEPF.NDEN 7 LAO AN6L6 . MA<IKC EUPs. At KNCte WHICH 

c ouadpan*»“ !T in 

cP-;LNP (L <^. LP .lO) = L ATANIUA ) 

IF (U4) 4i2.2E.;o 

?d ClJNTIN'JL 

TF (Ll> 2V,27,?7 
2^ :.P:LN^*(L^*LP .LOI sPI 

GO TO J7 
22 C'JNTlf.Ut 

I- (U2) 2a.24.24 

24 CPILNPILM.LO .L0)=rP3LNM(Ly.LP.LQ > ♦ PI 

2 a Ct NT iMJc 

CO TO 2 7 
2J CONflNUf 

IF < U.' 1 .£i ?. ?7 

eo -PjLNM^l^^.LPiLOI-EPSLNMILM.lp.LOI ♦ »I 


7 ClFiTlNL'E 

LP IE The tF^ECTIVC TIOAl LAG ANGLE IN DEGREES 
■ HP=Tn 7 L'NJM(LM,LO.L 0 )/r 

.iX GsQF = ( <L 1* »2 1^ ( 01 * «2 1 ♦ (E 2» *3 ) * <U2**2 1 )/U3 

SIjLUV is TML t TD. EWKOR UK YKMCLK.LP.LQ) 
a I oLUVst'.iOHV { a IGSOF ) 

s: 'iSQP= ( (c 2* *7 >« (U 1* *2 ) ♦ < E 1 •* 7 ) «( U2w *2 n/ (U3«* 2 1 
alGLAG IS TH- 5 TP. fcPOOR OK E^SLNM ( L« . LP .LO 1 IN CeCPEEf 
oI jLA*. = JjUKT ( f IGSOR 1 /F 


“C-FACrm “TSKCE CAFE TF"TFC AMFLTTUOr CT LJ/NITOLAC TirF 5 . S^E'^fMEEP. THE 

C LUNlwULAP Tia=.r. APr l.UMPEC* WITH ’•'Mr LUNAR '*'IOES IN SUiiPCUTTNr. CATA HENCE A 
C COPRlCTIQN MUGT EE APPLIED THE LUNAR TIDES WHEN T’hE LUMaULAR CASE CCMES 
C JH . SYMdOL ! CAlLY , ca:TDP =(1 ♦ SUN/MCCrj 1 «MOON s MCCN ♦ SUN. A« IT SHOULD FOP 
C THE LUNI SOLAR TIOLS. 

^ HALTOf -1 .OJO ♦ (LFNT5-U( LM.21) /( CFNTM.AI (LM.LP) 1 


TTCH? 


C YM^ALAO la THE LAG ANGLE OF TH£ 30LIC EARTH CONVERTED TC K RFQUENC Y-OEPENCENT 

C FORM 

YM=lm 

LALAG»OCOS< VM* Kt .*G 1 
SXLAG=LSlN(YK«tLAG) 

C aUbTRAcT OUT Tml SOLID EARTH ^inE 

SFAXt=TTI - XK»CTUAC 

IF (LF .EO. ? .AND. LO .EO. 2) SEA1=UI - F A C TOR * XK »C X L A G 
at Afi=U2 - XK-SXLAG 

IF ( LP *£0. 2 .AND* LQ .to* ?) SEA2=U2 - F ACTOR • XK ♦ £XL A G 

at A 3 = SrAl «*»? ♦ ?EA 2«*2 

SdA5=DSURT(«£A3) 

C KIND COCLAN. THE *'CfcAN AMPLITUDE 

CnC£Ar^=< 320. 037SS0C) •SE AbABl (LK) *Al (L-‘.LP»/ U .000 ♦ XKP) 

■ 'alTraOniHFl* •‘clY fSEAI* «Z> A I CS »42 > * f S 5 A 3 * v 2 ) T/SE A 3 
5>IG=CSOHT{ SI GEQR 1 

t SlGL la TMk 3TD. ERROR OF COCEAN 
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zz 

01 1 S 
0116 


ISN 

I SN 

0117 

one 

C : 

1£N 
1 .SN 
ISN 

0119 
012C 
0 121 
A 1 ^ 3 

4 C 

I SN 
ISN 

0 12 3 
0124 

49 

TTK- 
ISN 
ISN 
I SN 
TSN 
* Sn 

-TTzr~ 
0126 
0127 
0126 
0 126 
G 1 jO 

51 

44 

4a 

4 3 

ISN 
I SN 

0131 
0 1 i? 

45 

TTN 
I EM 

Dr37' 

0134 

4T 

C 

I SN 
1 5n 

t Clu 

01 M 

013' 

ol 

ISN 

0133 

5 0 


I£N 0140 
12N 0141 
15N 01h 2 
I £N 0143 
iSN 014« 
ICN 014c 
IbN 0146 


oi**e 

0|4»9 

I 5N 0160 
rFN’ Jlbl 


I iN 
IS^ 


iO 

'1»1 


^2 
6 3 


51 OC =( 3i.0.0376600) •SIG4RH LM)»DAeS (41 < LM,LP)>/( I *000 ♦ XKO» 

3IGSOf'S( ( 6 2» «£ )• (bcAl*«2) ♦ (El^«2 )«( 6EA244 21 )/< SifcA3** 2 J 

SIGPHb It ▼Mt. STO. EPRCQ OF THg PHASE ANGLE |K DEGWEE* 

51 GPMc = t)5QRT ( 3IGSON ) /F 

U4a»aE A2/5LA I ... 

FTWO’ THC^HASE ANGlT' OF ThE fJCEAN' T IDE . MAKING EUPE ME KKOA MUCH CUACOANT 
IT lb IN 

°HA5c=0ATAN( L4 4 ) 

.F (U44l 42«4e*47 

CPNTINOF 

(SfeAl » 49 » 4 7*47 

PHASEsPl 
GO TO 47 

•"■•■CCNTIKiUe 

IF (SEA2I 45*44*44 
PMAiEsPMAaE ♦ PI 
CONT iNUt 
GU TO 4 7 
CONT INUC 

IF ( bl A2) 45*47,47 
PHA5ti=PMA5£ ♦ PI 

•cD^mKlf^ 

PHAifis-PHAEt: ♦ (YM - 1 .ODOM (PI/2,0DO» 

injUkl That coolan is a positive amplitude 

IF (CnCFAri) 61*5C,S0 
CUC£AN=-C0C2 AN 
phaSE*pmA3£ - P! 

CUNT P,ME. 

CoNVcPT PHA5t FFCV DACIANS TO DEGREES 
PHA3E=OHASc/F 
IF (PHAiic) 'yC* U* 

PhA j£=t*MASEf 26 3,0 
CUNT IKUL 

IF (PHAEtll '>2, >5*63 
PHA5E*PPA 5£T 360*0 
CUNT INUt 

IF ( LP *LU* 2 .AND* LO .FO* 2) <"C TT 6^ 

PhINT jrrr L.M *p ,0 , constituent, DISTURSINC. body* Ft.LUUrNCY CEPtNOENT LOVE 
NUMtlct AI.O ITS LrtFOR* LA.* ancle AND ITS FPRCP* OCEAN AWFL’TUUE AND ITS 
PHAtc angle AND ITS t-RRUR 

Mf^ITE (5*33) LLL.LV,lPP*LOO*TID( ( L M , LP * L Q 1 • YK m ( l M • L F , L G ) . 6 t GL C V * EP 
i IC-LAG*C'JCL AN, 'ilGC.PMASE* SIGPHc 

I format (^IE*2a,A 4*4A. *M00N* ,5X ,F 1 ? ,3*2X,F10.6*2X.F<J,t,l**F«.e, 

* 4X,FlO,5,4:X,FiO,5*.''X,F5*4*2X,FQ 4) 

GO TO F-3 

' ■ XLVHUM=YKM(LM,LP*LQ} /FACTOP 


! SN 

0152 


61 GLOV=SI GLJV/F ACTCR 

I 5N 

:i£3 


4KITL 1 1 • 6fc ) LLLiLV, LPP*LGQ*TTDE (LM*LP *LU) * X LV NU M , S 1 CLC V * LF * S X GL AG 




1 fCOCfcAN.Si GC*nHASE.“ IGPH5 

1 SN 

0 1 G«* 

5? 

format (•>12*2X*A4*2X** MCCN4-6UN* ,3X,FlO,5,3x,F10*E*2X,P5.b*3X,F4,«* 




1 ^X.FlO, ‘j, 2X *F 1 0,6 * 2X*FO .4,2 X*FQ*4 1 

? EN 

0 1 B = 

5^ 

CONT IMUE 

I 5N 

0 130 


GU TO cO 

TSN" 

— 

62 

CXMTTF.UT 



C a 

ET TML Lf V- NUMOeP AND LAG ANGLE TC pQp JIDEa NCT tOLVuC FOF 

I5N 

015? 


YNM( LM* LO*LU ) = 0, ODO 

I SN 

01*- 9 


fcPSLNM ( lM, LP ,L )) =0, OCO 

I SN 

C 1 oO 


IF ILP ,rc), i .AND, LO ,tC, 2) GC TO ■» 

I 5N 

0 


-»«I7e: (cal lLL , LM ,lPP • LOO , T I Oc ( L^ , lP * LO » 

! SN 

Ole 3 

M 

FORMAT (•*i2*2x.A4«4X** MOON «,1?X,»,,,*»4<,»,,,*,8»,*,,**,^X, 




1 i:x,*««,*,>X,«,,.*,6X,*...**4y,«...M 

"Tsn- 

■01T5 - 


SO Tp 6 

ISN 

016 6 

7 

wF iVE ( » LLL ,LM .LPP.LOO.TI CE ( LM ,LP,LO) 

I £N 

C 166 

a 

forma* (4t2,2X•A4,2X,•MCC^♦^U^:•^l0X,•,.,•,9X,•••.’•SX••.•,••<3X* 




1 •,,. •*11X, •••*YX, ax, •.,.«,6X,*.,.»1 

I fN 

0!o7 

♦5 

CONT IMtJt 

ISN 

0 loij 

•>0 

CUNT INUL 

I SN 

Clr>5 


wW !TF (6,-6) 

1 SN 

0 f'C 

1 5 

FORMAT (lOX,////) 



C 

C # 

rt. Ah T * Jw r * • « • * 



C P'JT NA^’t 'Jf -IOAL CON6T IT'JuNT IN APCAY *»'IOE 

1 £N 

0171 


J'J 66 LM = I 

I SN 

0172 


DO 6 6 LP»l * 1 

! 3N 

C 173 


00 t'S LO» i • 6 

I £N 

017.* 

^5 

T ITL (L M,L»» ,L C) sKLANK 

tttt- 

~D 17 5 


TIDc(;*l*2)=T-lI2 

I*1N 

0 176 


TIDSC 1*2,3 )=T£122 

I£N 

51''7 


T:C2(^*i.i)=TS212 
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I 



w 192 
0194 


019& 
!3N 0197 


'!SN 

I 3N wZOl 


IhN 0202 
J SN C 20 2 


I 3N 0 20«» 


IfN 

I5N 

I 

I ?N 
I SN 

i sn 

I 3N 
■rsfr 
I ?N 
iSN 
ISN 

I «“N 
is«4 


020t 
02^7 
05C ? 

0209 

0210 
Oil 1 
Oil 2 

o-“r2 

0214 

021 r 

021 t 

Oil 7 
021^ 


TICE<2.2t£l=TS22£ 

DG 65 LMsl«2 

OU 65 LPsl • J 
YL»» = LH - 1 
LPr>=LP -~i ■ 

00 Ot L0=l#3 
LCO>LO - 4i 

IF (LF «rO. 2 •and* 

1 = I ♦ I 


LQ .60. 21 CC TO 56 


ISN 0219 

T^ N ' XTZ^~ 


15N 0-2 1 
ISN 02?2 


Xi>N 0 22 3 


U2*- VV< ! ) 

C IF ONLY TMl ilFFECTIVE L CVE NUMfiEP IS SOLVED FCRi PUT Tht EFFECTIVE LAC ANCLF 
C AND iTl, LRPJP equal TC 2L RO 

IF (NN(i> ♦ec» Cl U2»0«000 

iF (km;i .fq. c) E230»ooo 
C IF This TILc was NCT solved for, go TC 66 

TF tNNri3) ,E0, 0 #ANO, NNfl) ,EC. 0) GO TC 66 

U3 = Ul« *J ♦ 'U* ^2 

C FIND THU EFFECTIVE FrttGUtNCY-DEPENDENT LJVE NUMCER 
YM (LW, w*',LQ ) = 3SQ«T( U1 ) 

U*,=l»2/l 1 

C FIND THE f'FFcCTiVE F WEOUE NC Y- *) fcPENOEN 7 LAG ANGLE, MAKING SU^E WE KNOW WHICH 
C quadrant it IS IN 

fcPuLNo(LM,LP ,L0l=O ATAMU4) 

— IF (U41 32i3£,7T 
Jd CONTlN«iL 

IP iui 1 41 ,J?, 37 
41 EP-lNMLM, UP ,L1J =PI 

oO ro 37 

32 CGNTlNUt 

U ( UI ) 35, JA, TA 

34 LP.-“LNC(LM,Lf*,LO)sePSLNS(LM,LP.LO) ♦ or 

35 CONTIM*- 
Oj to 37 

33 CONTiNUL 

IF (U2) 30,37,37 

Jo LPSLNS(LM,LP ,L0l*rPELN5<LM,LP,L0) ♦ PI 
J7 CUNT iNUti 

C tP 13 7Mf LFFECTIVE T 1 0 *L LAG ANGLE In DEGREES 

LP=EPiLN3<LV,LP,LO)/F 

KiGi:or=( *2) ♦ 

C 3loLOV 1^ ThF ETC), ERROR OF VKS ( L M , LO , L 0 ) 

5IGL0V-P5nHT(SIGSQRI 

•^I GSQf =< ( fc i*’ *? >' <U 1**2 ) ♦ ( F l*«2 I Y ( IJ2* ••2 I I / (U3*w 2 I 
C jiGLAG IS lot STD, C.RMUR OF EPSLNS ( LM , lP , LO ) IN CEGKCtS 
SiGLAG»C*SOHT(i£GSQRI/F 
C 

C DO TmF. CCLAM ’’’ICLS 

C “ 


ISN 02?4 


ISN 022S 
ISN 0226 

tSN 02 c7 

“TSN orre 

I £N 0219 
ISN 0230 

ISN 0 23 1 
I;»N 0 2 32 
ISN 0233 


T«N or34 - 
ISN 0I3F 

ISN 0236 
TSN C^’7 


TSN 023F 
“TTN 

T PN 0240 
IfN 0241 


c YM^KLAi u The L AG ANGLE OF THE SCLIC EARTH CONVERTED TL f- cEGUENCY-CEPENDENT 
C Form 

CXLAG^DC 13< VM*XL AG ) 
iXLAG=r>i'»I(4( Y«^*XL AG I 
C SJtJTRACT CUT THE ^Ol I D EARTH TIDE 
6r:Al=tJl - XK«CXLAG 

- XKW'SXLAG 

SEA25s‘ t A I « •* 2 ♦ StA2**I 
Sfc A3=DSCR^(St A3) 

L FIND CuCEAN, TH£ nCEAN AMPLITUDE 

CUCl AN=( 1 46. 77 ;oOC I • SF A5* 31 ( LM ) *P ( LM ,LP) /( 1 .COO ♦ YkP) 

GaiV* =< ( r I ^ *2 I* I Sn Al* 42 1 ( E 2» *2 I A ( A 2* Ag I ) / S fcA 3 

M ,= OIOKT( GSQR » 

C j lie IS 1 Mi. .JTD. i RKoR CF rOCFAN 

ST GCst TAG •‘^T95nTn*SIG ?Bl TLM) ^ABSTP(LM,LP| I /n .OUT) ♦ »KP1 
iIGL0r=((E2*'*2)'(ScAl*42) ♦ CE!««2 )«<SEa2a»2M/(SEA2»*'?) 

L olGPHS IS the STD. ERROR OF ThE PHASE ANGLE IN DEGREES 
»I GPHi=D5 iH T ( STGSQR ) /F 


U44 2 3C A 2/3fc A 1 

C find THL PHASL Ar4.",Lfc OF THE OCFAN flOE, MAKING SURE Wfc KNOW WHICH GUACRANT 
C I T IS IN 

RHAi£»DATAM U44) 

— fV^Al S3, *5,51 

Dd CCNTiNUt 

IF ( 5E Al ) '3R,£7,57 
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1 


1 


1 


I 

I FN 
I *:n 
I FN 
I tt- 

“ ^ * 
0 'A * 

0 0 

02CI- 

39 

o £ 
n«* 

pma jc*PZ 
G 1 r ) ^7 

CtiNT iNUf 

;F ( OlA • I -3' * :'*• ‘'A 
pha : t=pMA jr ♦ PI 

TTN 

0?f T 

*55 

CciN T INtJt 

T SN 

•02WC 


O'J TO *■- 

I SN 

C 2A9 

bO 

CL N7 

I ^N 

r. “* » r. 


l^ IJt. A.’ ) ft) . i 7 

I FN 

OrFl 

5a 

PHA‘^£: = P^Ari ♦ ‘‘I 

T ^N 

02-’£ 

^.7 

LUNT INMf. 

IbN 

0 ?'. i 


PM A »F.s-*N‘A wi * ( YM 



C iNa'JPi. ThA*' cDClAN I*. 

T^'rj 

C7’"3 



i ON 

C J' 0 

73 

Ci f.LAr;=-c JCf an 

T'.N 

0 i •• 


PMA»L=PMAaL - P* 

! FN 

0^07 

74 

C( NT :^ ML 



C CLNVlPI P*-*AS^ FLf-^ kA 

I FN 



> £sPmA 3c/F 

1 a^- 



ir (PHA.'L) 7c*7c7 

I Orj 

C 0 

7a 

PHA:'. L-PHA jL ♦ 1- 0 • c 

ITN 

orr I 

7f- 

CjNTIIP»C 

I b^. 



ih 1 pMAr .) 7 7 , -»t , » 

I FN 

• 

77 

PriA c. =PF A br ♦ It - •’0 

T c». 

u 

7« 

cr nv it.^L 


•• ■ r i’,! iT ..'in Lii.p.nt ctiN^'; TUKN’’, or -.TunnrN". mupv, e = l-.«-_kdent lgv‘^ 

C A* 0 l“t» LAi ANOLf AN*' *_«W0*^* OCE.AK A » TU'Jfc ANO IT$ 

C J» • f*HA‘c ANuLt- AN.') !Ti, rp:^rp 

ISf. 02^- *-•*"'£, («•'!) ILL t5_ *»# Lrn,i OT, TI CcI CL« .t ’ *L0 I .VK\>< LW, Lf .1 <. ) t ■»! *L0V . EP 

' - i,!i^r,LAr.*C3CrAK.^I:.<',Pr*A3r*-!Gi'H^ , ^ c 

I?N O?*-'- 71 r)><MAT (Al2,.iX,AA,;x,».SLN»#o***^lC**t*2’<*Pl««t3»2**^»» • ■•) **5* 

I 5N 0 2- 7 us) 1 Ti t •. 

,N ^ THr LL Vfc AND LAG ANfiLf ''O 7S£Pr rDfJ NCT ■(<LV':0 FPP 

^ Zt-. GZ*'^ ’ YKKL^'.L'^iL'J)* :• rOD 

ifN 027C rPtLNV (LM, L^» .LO ) =0. ODU 

ItN C2**l I'^E C't=^3) LLL .LM ,LOP»Lr.C .TTCr (LM ,lP#LQ ) . ■ 

fCKj 0e^3 PoPMA" (4l^*2X*A**i'>*#*jUN**\ E*#******^^*******^’****** * ^ • ••* • 

1 I ; X »••••••' X »•••••• 6 X f ••••• #bx t *•••• ) 

O^?"* CD CJNT iNiit. 

: :-N C27A <f,0A) (tatcj)# 

I^'N 0«7C 4A ^V>LMA‘»* (////✓// *£0 X t 1 A*5 * A -* ) 

C PSilNT < L^ Lf V(L)# ANJ XINCL(L) 

I TN oarf. DO I 0 L= I • »o 

TT=T(L) 

! TsN 327A CALL AP-.>(TT) 

1 FN 02*^G jJMsC.OO) 

T SN C'PC 7 0 L'4=l*J 

!SN 0261 YLO=u- 

tSN 0£>^2 r>0 7 0 lp»i»3 

I SN 02*^2 YLP = Lf - i 

^ xrnrrUL l.LP^LO ) =VK ‘Jif L «• LP *L 0 1 • Cr^.Ts^ • ni ( U M> * f (L A , 2 » * A I (1. , lP> *CL PQm 

I ( i^u , LO » • VL V * Ahu'J*' ( L *1 ,LP • LO ) / ( ARG C' P ( L *' • LP • L Q > “S S ) 

1 cf^ 026c Ap'jNU'^.LP*^ C)=VKS(L'x,LP*LO) •CrNT‘^*Hl ( L 6 ) t C ( L . 2 I - C ( L ■» • LP ) • GL P Go ( L 

lP*LO)'YL?^*‘AhOC:(LM«LPtLO)/(APGDT''<L»»tLP*LO)*SG) 

I 5N 0267 5i-»s:>l.v ♦ X‘UJCK(LM,lP.LOJ ♦ XT LN ( L ^ *L 0 ) 

I SK 70 ClNTJMM 

»‘“.K 0£*aR XiNCL(L)='.rU*6 ♦ xaD 

! SN O’-JO lO CUNTlni).. 

rr?i wwi’x (w.i'5) 

I SN 030.- ITc ( «• 1 00 ) 

I S.N 0293 100 ru~f*AT (V3X*'TIDAL I KL t N A TI cr- • # // » 

I«N0 29A <iRIlf('*iO*) 

!SN GZ'ft lOi hU^-^’Al (i5Y,*TI'*fc*.|6Xi*I NCL !N A’^ ICK* ) 

J aN! 029c *Kli'L(^.»*»02> 

IPN 0297 lOc (jJX*«(JAN 1* 0 MRS U7 » 1 ) • . 1 2 X # • ( A PC b t C > * # / / I 

J tU Z 2^ e DO 1 2 L t ^*T 

*T5rJ 02-7T aRTTF: T^.17) L .T CLI • XT NCL (I- ) 

l<zK 03J0 12 CONTlNUt ^ ^ 

ISN 0301 IJ ^JP'-IAT < jC X t • bX * lO * 3* 5 • 1 C X • Dl 5 • > 

C PPINT RESjlO'JALi 
I FN 0 30 ? *R I rt ( *5. 1 9) 

I SN C 40 3 <•» ITE ( t'* I 03 > 

I 5N O30A 100 FU9MA1 ( A6X, •Rf.S IDUALS* • // 1 

tin ^33? lOA FORM AT * { 13xJ *T I«E* . J *3X, • PE'jIDUAL • ) 

!SN 0 307 aPITL {'>.1021 

1«:N 02C2 aUMilG = C*0 






; SN 

0 » ^ s 

J«J 4 J * 1 t ML X 

I J>H 

- n c 

rT=T^xP( J 1 

I «N 

0 3. 1 

C ALL AF r,( tt ) 

I SN 

0 •*: £ 

^ViM= “ , 

T 

J *i ? 

3 k ,i 

I SN 

331 4 

VL MsLM 

I SN 

031 f 

f^i 3 LP- 1 • ? 

:«5N 

3-1 * 

VLP=Lf - 1 

I SN 

Oil 7 

DU ,• 1 (J = i , J 

T ‘-'N 

3 

AMU 1N( L ^.LP. 



1 (L» • LO) • YLM. 

I SN 

C ’i Q 

(LM, ,i_ 



JPfLO) VL-i.Ar, 

TTN 

C 3T^C 

♦ xM 

! :sN 

0 1 :. 1 

J C-NriN»Jf 

1 SN 

C **1^ 

3UMs:,u‘> ♦■ XI 

I SN 

: 3£“ 

PLi IsxlLX M J 

I SN 

03^4 

SOM ^ lO^t-c; i. 1 • 

I SN 

0 *■ 

ITfc 1 o» 1 Jl 

I CN 

6i3<» 

4 CuM-inUf 

T SN 

C 3C 7 

XNsMtVp - : 


-V 17 S- 

ifN 

l•^ 05JC 


» SN V i ■« 1 

ISN 02 I Z 
1 SN C 5 J J 


• UPfLO) •C*=^NTM*dl (L'4)«C ' *\ (L***Ln *OLPQM 


• <L*i',LP»LO)»'CP' 4 T.">»^l(L*<)*C{LM, 2 )»F(L*'« 
*.Lt,La)/<APSCTe(Lw*LP*LO»*«S) 

1*l-P#4_0» ♦ X?UN (L«#LP*LQ) 


LP i •’■*LPC5( L 


J • Tx. XP< J) i*t «T 

-iirMAs^nPTf siiM^r '/XK> 
i*K iTc ( o* * c'l ) A 

^ ^ (////✓, :^ «Li -,MA = sacT((‘LM jr cauAKc^*/c p£xf-i »» = •, 

“LOT C:UT XiNCllL) Vi. T(L>* nFSirtCAL* CN COMPUTCK PAtifF 
CALL PLf TcP( MT .Mc«P> 

''L'Uri 

cNi") 


56 


I 


I 


1 


I SN OOOi iOHP'JOTINfc «»L1JTE W ( XT . ME « P ( 

C 

c 

C THir. iUnROlITiNE PlOTi » CCNTINDOEi C'JRVE PLUS E«PEe!XtNT*L !J»T» POINTS 

C Vi. time tie), it also PLOTS THf ^esIOUALE VS. TIME T(ul. 

C Time INCfiEAic.3 00*A TME PAGE AND TPE KlDTp CF THE FLCT Ig lOI SPACES ACECSS. 
C theme ale MT points in The C'JRVE. WHICH IS STOPF.D IN AEBAT IINCLILI. ANO 
C Mr.XP lXPLFIMENTAL SA’’A POINTS, S'^CEET IN APPAY XIEXU(J). THE POINTS OF The 
C CURVE APE printed AF 0*« ANO THT DATA POINTS AS X»S. 

C 

c 

I SN J JU_- ■ “IPPl rCIT PE AL t s I A-n. E- J ) 

I EN O iSA JIMFNiUlN EESIOliJ)) 

liN :J3E OIMeN.ION BINrltOll 

ITN CJ3C OIMf.NilUI TC XPI ESO I , M eXPI ?00) , t (7 00 ) . XI NCL I TOO » .NT CACKI TOO ) 

I SN eOCT CUHMON/HLXC/ T£ XP ,X lEXP.T, X’NCC .NTPACK 

I SN OOOe CAIA CEE.EAf .HLANK/l FO ,1HX . 1 h / 


C DeIEBHINE the width HP THE PLOT. HY FINDING THE MAXIMUM ANO VIMMUM POINTS CF 
C THE CURVE AND the CATA PCINTS 


rr;N 

'7ZZ^ 


XMAK l=XlNCL( 1 > 



r «N 

O'jIO 


XM IN isx:ncl( 1 ) 



isr4 

1 


OU jO I=2.MT 



: SN 

0012 


OsXlNCH I » 



I .^'N 

301 


IF «XWAX1 - P) 22«21»<.! 



! CK 

00 1 * 


KMAX t *P 



I SN 

3 JJ S 


on ro ?o 



1 SN 

0 C 1 r 

> > 

CONT iNUr. 



T5N TJi T 


xi.(P-xMTNn2A,.Tr*33 



T £N 

001 e 

s, 4 

XMl\i=P 



ISN 


JJ 

CUNT INLL 



!SN 

0C30 


IF (M-xP .r.j, OC TU 4 



1 dN 

0 32 P 


GC VO i 



! SN 



XMAXW'MAXI 



I SN 

0 3i4 


XMIN=XMINI 



I SN 

OCIf 


GU Tn so 




' " o 

ciMTiNur. 



TSN 

332? 


xmax2=«ICxF( 1 » 



* iN 

3 32e 


XMiN2=Xit XP( 1 ) 



: SN 

: 3?c 


CiJ 4 3 J-i.^cXM 



! TN 

0 o?c 


r>s X I£X»J{ j ) 



iSN 

30^1 


IF IXMAXi - P) 42,41 •«! 



I SN 

C 3 


XMAX 2 = p 



I SN 

d:**" 


GU *’•3 40 



TTK" 


' i,"! 

CCN^INUt 



2 SN 

3 325 


IF (P - XM.‘N2I 44*40*40 



: JN 

00 

<•4 

XMlN2=P 



! f.N 

CO. *7 


CJNTir.MJc 



2 SN 

03T-? 


|F (XPAX2 - XMAXll 4£»46,4n 



I SN 

OC 5C 

40 

XMAX *KMAX2 



I SN 

0 0«*0 


,0 TO 4 7 



t 5N 

0 “ i ; 

4 S) 

XMAX SAN A XI 



TrN 

3 0^ r 

♦ 7 

C^NTI^.Ur^ 



I SM 

00**2 


IF IXVIN2 - X»*IM) 4e,44,44 



I3N 

:;-;o 


XM 1N=XMI .42 



I SN 

c S 


\j j 7 0 53 



I s^i 

0 c - * 


X V ( s n M I ,N A 



! sN 

3C«7 

•iO 

CUN7 INML 





C C04PIJTT FAC. jO that ALU OF Tne OCIST'; FI* 

rjN paplf 

I «N 

c « 


F4L-U0,J.j''/(xmax-xmtm) 



TSN TTC* ? 


^^rrc (6, 3) 



! SN 

OChC 

3 

rOPMAT < IM 1 > 



1 SN 

u3:;l 


SCALC«XMAX - X'ilN 





C 4:^ITe the lLW-IST hICmF^t PCTnt* ANO 

OlaTANCe 

C-tV'AFfK T*“cM 

I SN 

0A5*"* 


APITE ,13I» XMIN . XMAX, SC Al £ 





C 3=T 

APPAY initially eOUAL TO MLANY. 



T ^N 

; 3f: 2 


OfJ 26 J*i • 1 )1 



T zh 


?.U 

t1 1 Sf ( J ) = )L any 



TSN" 

3 z 


or 23 L = I ,MT 



I jN 



J4si 





C (^JT 

CUilVt POINT IN LilNElJJ) 



: ~u 

3 I'."' 


JJ=(x2NLi_t L) - XMlNMFAC A i.Si 



1 

)3- t 


ri4Ni{ J J )-J£E 





C 3r ?i 

IF AN =wP£PIM-n1AL OATA OCIn7 SHC'»L0 JC 

PLUTTeO, 

HY CPrCYINO TPP VALUE OF 



C NT'S 

AC^ ( L J 



I5N 

o C .• * 


If INTFACKIL) .cO, ?) GO TO 




c e 1 


NNsNTf At Y ( L ) 





C F JT 

r.XPCPI ml >JT AL 2A*A PCINT IN FINE<J4) 



T 

C 0 62 


j4 = ( x::.xp( NN » - XMINJiFAC ♦ I,5l 
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I 


I 


ISN 

l?.N 

I £N 
ISN 

00c3 

0064 

OOt « 
0 OCc 

t31MECJ4|s£K3 
60 CONT INUL 

L print the number of the curve POINT, *IINE, ANO THE N^»MPCR THE CURVE POINT 

WRITE (6*29) L*(UIN£Ol>« Jlsl,lOl)»L 
25 format ( 0X*1 S,2X*1 01 A1 ,1X, If. ) 

I«N 
I EN 
1 SN 
1 SN 

ISN 
I SN 

0067 

ccoa 

OOoO 

0C7C 

s7072 

CC7i 

C PUT 
29 

C PLOT 

APPTkV DTNE(J) P4 Ck TC Blank 
BINE( JJ IsULANk 
bXN .( J4I sUlANK 
LUNT INO;- 

IF (.Mt.xr .SO, 0) GC TO 130 
RE&IUUAL3 VS, TINE 
DU 132 J3l,MExP 
RtS tO( J la 7.0 

2 tN 

TrT4 ■ 

l J2~- 

CCNTiNOr 

I SN 

0C75 


OL 1 Oi Lal.MT 

Z SN 

0076 


IF (NTRACK(L) ,EQ, 0) GO Tn ICI 

I SN 

C076 


NN = r4TRACK(L) 

I'.N 

0079 


»:.3lC(NfUaXIEXR(NN) - XlNCL(L) 

I SN 

7vdO 

tot 

CuNT lN‘if 

1 SN 

0C61 


XM AX2*Re 0 ID( 1 ) 

I 3N 

0C4i 


XMlN2sRP‘ilO( 1) 

1 SN 

J C b J 


DO l4C“Js2,MExP 

1 5N 

0C64 


R=RE6I0( J ) 

I iN 

0CS5 


IF (XMAX2 - P) 142,141.141 

I SN 

OCMt 

142 

XMAX2*^P 

I SN 

OC 37 


GO TO lAO 

I SN 

OO-JtJ 

141 

COUT INUE 

1 SN 

0C69 


IF (P - XMIS2) 144,1 40, 140 

1 3N 

0070 


xmxn2=p 

TTR" 

JR 1 

1 40 

CONTINXjr ■ 

t SN 

00-->2 


oCALt»XMAX2 - XMIN2 

lEN 

0093 


ifeRiri. (6,3) 

I SN 

C07A 


• RITE (6, Mi) TMi n2. XMAX2 , =;C4L E 

I SN 

0 c ? 

1 J1 

FORM A*, (/////,10X,*RlNlMUMs*,ri0,o,l0X,*RAXlMUMa*,FlC,S,l0X, 

T SN 

00 ;e 

1 

• ECAL C** ,F 1 C, 5, ///////) 
FAC=1OC,0DC/(*max2 - XMIN2I 

Z 5N 

CO^ *» 


DJ l2r Jal Cl 

rrrr 

•TOTS 

IZ-j 

JINC( J ) =I1L*NK 

I iN 

0079 


OL 12R Lsl,MT 

I SN 

0170 


J4=i 

! SN 

0 17 1 


j j =l -XMI Ni )*F AC ♦ 2,51 

I SN 

0104- 


riiNc. ( JJ ) =uE£ 

I TN 

0102 


If (NT«ACK(L) ,E0, 0) GU To i^^C 

lEN 

CIO® 


NN-nTF ACn (L ) 

I ‘‘F 

0 1 Or. 


J4 = (Rl 5iU(NN ) - XMIN2)«FAC ♦ l.Sl 

TEV 

• OlOT" 


TTr'irrjiia&Kr - 

I SN 

j IC ® 

1 oO 

Cunt inuc 

I SN 

Oi: 9 


• (L,2S) L,(fHNE(Jl), J*al,lCl),L 

: EN 

oil 0 


RiNL( JJ)aMLANK 

I SN 

0 21 1 


dINr.( J4 )sHLANK 

1 SN 

2 ?N 




on: 

I 30 

CONTIr.Ot 

I SN 

> n A 


R£ "URN 

“TEN* 

Til'S 


zup ■ ■ “ 


C 

C 

■C' Tk<I 3 TJJMttklfJTlN* ^“AC.3 In and "^PTnII OUT The satellite fLS’^ENTS AND NODF 

C »aT£, all INP'JT U. OEGRFFa ( GK< DEGREES »EP DAY) 13 iw^fOtATELV CCNVERTEO ’’O 
C PAuiAND (U>) HACXAN3 PER DAY) FOR LEE IN ThE FPCGRAY* 

C 

c 

implicit me AL«-» ( A-m, C-^ ) 

DiMi. Nt ION bAT ( 1 A ) 

COMMON/^LKD/ >«LAv.. 0 .r, MtGAM, OLLMN .CMEG A3,OLLrUN*ES #*!3.A, :.XI. XNCOEf 
■" 1 'TZCirr* VI M30N • SON * TmS , xn^m ,x ncs* F, /meanm , X wF ane, X , xM.ioTS • 

2 OMOJijT »MEXP 
COMMCN/PLK 1/ S AT 

R=A3 (b.T» (3AT(j), jsl^lA) 

7 FORMAT ( i 3A6t A’ ) 
wR I^fe ( n» 2 > 

2 r JRf«Al t iHl ) 

**hITE <Otrt) (>AT(J)f Jsl.l4> 

S format (/////. 1 ox • •satellite* . I X ?A'i*A^ ) 

READ (j.n A,ttOl,Xll.TZ£FSt<NCOEl 
I Format ( JFIC .£) 

xi=t ii*F 
OsfJl-*F 

XNoOc* XNODcl *F 
wP*Te (c,J» 

\ format <////// *^0X •• SATcLL* TE data* •/////» 

- "d» ITt ( 'j,'*) 

A FORMAT C 2AX • • A * , I A X ♦ *E • • I cx • • I NCH NAT IGN* . AX . • Nl'Cc n AT - f , x , • T Z EPQ 
A*tllXt* N JS c * ■ / ) 

WRITE (oi3) 

•3 RQFM AT (?.'X**UO**E CM)•,21X*•<DEGCE£D)••5X.•(0EG/CAY)•,‘^«♦•(CAySJ 

I • • “X t • < OkGPE F3 » • •/ / > 

wA I Tt ( Nt o ) A ,xl 1 # 01 # TZERC. XNODE I 

N FORMAT ( 2yX, N(F 1 3.A, 5X ) I 

iTETUWf" 

END 


I SN 0 70 2 


ISN 0003 
I EM CwDA 
1 SN 000£ 


:tn jooc 
ISN 700> 
i $N :co® 

1 <N 0009 
: SN 001 T 
I «N 001 1 

TTN 0712 

1 3N C )3 J 
I £N O'O: A 

: EN o:i E 

I 5N 0 0 1 ^ 

I SN 0 71 7 
! fN 0 312 
ISN v71R 

TEV'TOTT 

I'N OOZl 

: f N 0 0 2 2 
I SN OC 3 

ISN 0 72-* 

I SN ZCZt 
-TrT*rT76~ — 
ISN COOT 
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I 


( 




»3N 0003 ^UH**'"uTInl 

C 

c 

C TniS KcACS IN ANO nWlNT? OUT Tnfe PC^^ITirs ANC AK'‘fUL AO SPeCC* OA 

C THi 4UUN ANi •‘UN. INPJT *N OEGWfrf (C-* ’5C'j»^Ltt »»t0 J AV ) |‘ MMEDIATCLV 

t CONVcR^IiJ TJ kAOXANj (OK KACIAN'- P£K DAY) FfJO USE |K TFE PRUGKAM. 

C 

c 

f ;N 0303 I*A»'L|Cn AL»-i(A-K,C-Z) 

I IN 00 3 4 CummCN/OLKH/ «L AG .Q •C'LOAM, OeLMK* CMC'l as* 1 EL 5 UN.Sf. •KI 0 *A«F#XI. KNODE • 

A TZ::'5'’fXl MC'UN, X 1 iUN •TM0,XNnM,yNO«*®**A'EAKM,K*^fcANt*TM'0(. r«^» »MOCT« » 

■ ■ 2 (jMr,DCT,‘Atxc» 

IfN :rOt <?.ll IMt. ;a 1 , Of lT Al lO'AgGAa t O-L*^ A2 , XNC »»i •C'^COT I 

! aK OtOc R£AJ X ID 1 • XI 02 *T M's* X^E ANt • X 'AC an2 

ISN C*3T I FUOMA*' (JF10#*I 

t X I M J^NsXt J I • F 

! $N O vO** X I s.MN*x 10 2«F 

!iN 3010 DtL4N = Ofcl.f A1 •f 

I 3N com i^rL iLN = JUL TA i -f 

co;? CMC‘;A»^=-:.At:^Ai r 

’£N 001* CML0A;;*1 tl£GA '*F 

I SN 5: I 4 XHCANM=XME AN 1 *F 

! ;N OClt XMLANE =XM ; AN2 • F 

I F.N OOlc XNJM»XNL)M1«F 

ISN COiT OM jL>OT*UMUCT l'*r 

I •:n co*p xn:je X- 0 . 4"» 

! SU 0 ox ? iTE (6 . 2 ) 

'TS*r'wT 2'0 d' rciPMAT ( 1 ;y, /////,« ^x, *mcOn ano tun data* ,////! 

I5N 0 0? I «K ITE (F • ) 

13N 0C2f 0 FuRMAT ( J J » . *C MC'iA • . I 0 X • •DCLL TA • . 7X , • I MCL INAT X CN* • 7x . • TX'i ••ax« 

1 •McAF ANO^'ALV • *JX , • Mt. AN MOT | CN • , 5 X * • NOOE FATE*./) 

I«:n 002 2 «K ITC ( o. 4 » 

I £N OOcA *. format (^IXiMi. G^Fs 3 ) • • oX.* (l'£Goe = s ) • .C.X. • ( CFG»^tt S ) • t 7 X, • ( c Ars ) • • 

1 * X . • ( DL«j»^EE i ) • AX • • ( OLO/CAy ) • X. • ( f>F*j/i; AY » • .// J 

ibN 0 02 5 aPITF (t.,0) C»'~0\l»O6LTAl.xi'01.Ty0*XMCANl*xNOi*l«CMrCT; 

Tt:t 4 C 22 P r. FHFMaT ( oV . *MnrN» , t ox, 7 (FX 4 ,ex ) , /) 

!SN 0327 (*,♦») C'*‘^f.Aj,U\;i.TA 2 ,xl‘ 0 n .TMS.XMfANi, xnDS* 

ITN u02fc o FjkmaT (ox. • S t J * ,; :x otFlO.A ,bX) ) 

I«-N 0 32F RETUWN 

! SN 0030 ENP 
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oor»nnnr»r|rw>r.r'nnrk 




>lMi40tir INC ht Gat- S( wT ) 


Tf«T?l CAO^ 

Th; i . 'JH <0UT IN ■ la 4 
PA\*t3 4D» - 4 07 JF THt. ! 
VFNalafi 111* f»«n jK 4MMF a? 
(4Js»U‘»T 

Vmc. p^.*ipuitf Hf***-. I t> 
aLiatasICN* (2> PtAO j'iH 
C4 LC.JLATc M:.ANS* ‘s^'ANDAW 

crr^iciEN^r* •^-values * 

(4) ai, iNT TM£ cfc^iiLTS* 

rJK jCurutTiNr. lc 
TME NHMt^Ca CK IUjSFQV 
Ga£A’'£r. Than m**, «me 
!*• JaLUitT SfcLFrTIUN 
MuuTiPLC rruR*i‘i!(N* 

AFTfP N^rjCNlNG TRCM 
rrSTCD 70 CHcC< a^-cthep 
CiJM^Awen AiAlNFT A b*<ALL 

singularity, 

TmIj ^.UiRJUTIN** call 

.TA»'A>* OnCE^'f *^iNV, MUL*^ 

APE aTANCAPU S UHRrNT I NE S 


IE*: nf •'HE Miit^IPLE PE'jPES'iXCK ANALV3!^. 

MOOIFICA^’I IN IF Tnr PRCGOAM L«PLA|N*iC ANO LIFTED CN 
CM ^v«rEM/7oC SCIENTIFIC S»JUR01’TI\E PACKAGE, 

MANUAL* PPOCRAA NU'^REW ?oCA-CA-0?A* FiFTH ECITION 

TH <iJ WfcA3 The POOriLtM PARAMF1E6 CALO FCP A PLLT|PLE 
SET ^ELECTXJN CARi*)", (3) CALL Tnt 5U600UTIKE5 TO 
0 DEVlATIQAS, 5I»»f*L^ AND pULTIFLE CCOUgLATICN 
thn ANALYOT*^ nr VAFTANCE rr»R MULTIFlE OEGOESSTCN, 

NO I5> TAVC Tm= PLOOESSITN CCLFFICIfmS ANC TfCIP 
VNUM, 

ANIONS, Mr*P, IS CALLED K IN ^Hli dUt^ C 'JljT I NE . N ML£T 

Ff M IS ▼HE njmijER of variable*:. 

Capos ape not opesent. The procnaa <akn#)T perform 

SUORCutINE mINV, The VALUE CF Cc^LR^'INANT <0ET» |« 
THtt CC^RELATICN MA'^PIV IS STKGL'LAP. SET IS 

CCNSTANT, This TEST MAY ALSC FE L 5 E 3 To CFECK NEAR- 

S SL’OROU*!NC£ CO^kt (AHiCh, IN TU^N, CAlLs SUfcROUTINF 
C, ONLY SA'^A !> FOUND LISTED H*«t, SINCE ThE OTHERS 
CCnTAINEO IK THE machine. 


" THE SLlTTiT aTLiCTirM CAROfSI GIVES <• ThC MJmBEC OF CC^FFicIEnTS SOLVED 
FCh* AND (I.AVC(J), WHICH CJNTAINS THE TIDAL CCMiITUKNTS SCLVED FOP* 

THc table rttLJW GIVES THt TIDAL CDNSTITJENTS ANJ Tm£ 1 P ISAVtlJ) VALUES, 


L 

M 

p o 

f I CF 

01 »,HuDy 

1*' 

• AVE(J) ' 


1 

0-1 


MUON 

1 

AND 

2 

2 

1 

DC 

ni 

MOCN 

3 

AND 

4 

£ 

1 

T 1 


MOON 

« 

AND 

< 


A 

1-1 


MOON 

7 

AND 

e 

i 

1 

1 0 

K1 

M0CN4 SUN 

9 

AND 

1 0 

£ 

1 

1 I 


M'JCN 

1 I 

AND 

1 2 

< 

1 

i-4 


MOON 

13 

AND 

14 

2 

1 

i 0 


mjon 

15 

AND 

IK 


1 

e 1 


MCCN 

17 

AND 

1 b 

ZT 

^-1 


■ TTnCFT 

T9- 

AND 

'2V 

£ 

2 

0 0 

M2 

M 7CN 

21 

AND 

22 

£ 

£ 

Q 1 


mcCN 

23 

AND 

24 



1-; 


MCCN 

25 

AND 

26 

> 

£ 

1 0 

K2 

MUONaEUN 

27 

AND 

78 


2 

1 1 


MOCN 

29 

AND 

3C 



2-1 


MOON 

3 1 

AND 

3? 

* 

« 

2 0 


MOON 

"■ • 

33 

AND 

34 


£ 

2 1 



'JS 

ANP 3 6 

£ 

1 

5-1 


SLN 

37 

AND 

36 

2 , 

1 

3 C 

PI 

F'JN 

39 

ANC 

40 


i 

5 1 


SLN 

41 

AND 

42 


i 

1-1 


»LN 

43 

AND 

44 

i 

1 

1 0 

Kl S 

SLN 


• 


£ 

1 

1 1 


SLN 

45 

AND 

46 

■“Z 

T 

"2-1“ 


-^UN 

47 

AND 

4 3 


i 

£ 0 


SLN 

49 

AND 

50 

> 

1 

£ 1 


SLN 


AND 

52 

i 

2 

C-. 


SUN 

53 

AND 

54 


2 

0 0 

32 

3LN 

55 

AND 

56 

£ 

£ 

0 1 


SUN 

57 

AND 

58 


2 

1-* 


SLN 

59 

and 

50 

2 

£ 

1 0 

K2 S 

SLN 


- 





“ SV7T 

61 

AND 

6Z" 

£ 

z 

2-i 


SLfi 

53 

ANO 

64 

? 

£ 

2 C 


SLN 

65 

AND 

o6 

2 

£ 

2 1 


SLN 

67 

AND 

56 


C 

c 
c 
c 
c 
c 
c 

c 
c 
c 
c 
c 

C EACH CJNSTIT'JENT has Two ISAVF(J) VALUtS, SINCE «E SCIVE FQ« TWO 

c UAPAMdTERs A LOVE NUMBER AND A LAG ANGLE, IF <nE wISH TC SCL*"* FCP ONL Y A 
■C“XOVT''N*JMTJFP“ (ft3pcTNG“THF LAC ANGLr TC BF ZEROl , "THFK C^LV TFC TTP5T xSffVEIJT 
C value of That constituent should appear on the subset ScLuCTICN CARD, 
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IfN OC .03 


Z zH 000<» 
I3N 

ToN 3 U 0 e — 

liN 00 0 *^ 

i *jn coce 
I c cc 
X SN OCl 0 

Il^N 0*>1 1 

I 6N 0 0 I £ 


1 SN OUl J 
I FN 001 A 
I 5N C ' 1 c 

I OCl e 


i:n coi"» 

I :n 0019 
l«N CCl^ 

\e^*4 ooac 

1 «N ODi-l 
“TofT VT^.Z 
X 0 02J 
I£N ;02A 


XSN D 0 i£ 
T r\ ooi«i 
:?N 0 : 2 7 

IFN oc-?e 
ISN :c.’» 
I ?N JOJU 
TTTj nirrr 

XSN 00*2 
I SN c :• 1 j 
1 «N oc:a 
I^N 0Ci*5 
1 S^’ 0 - >c 


1 es 0 : j7 

Tcn OC.ld 
I f N C';!*; 
ttr; OCftC 

xeN 0 C 4 X 
I 2N 00*2 
I*“N OCA 4 
I C 0«*«» 


: 5 N 0 T-*e 

I £N C 0-*4 


I 


I 


C 

C THfe FULLOlilN'a 0 1 NS I ON t» NUST BE ^RSAT€R THAN CM fcO*>Al. TO T*-6 NONOEP 
C VARIABLES* •* 

C ___ 


OXME N< lUN XBA^> C 701 *b TD( 701 #D (TO) *P y( **0 » # ISAVE< 70 ) # V( 7; | , 

1 aii( 7v) *T( 70 ) •«( 70) 

C THt rLLLO%lNvi niMtNSION MUST BE 'OrEATER THAN OP LC^AL 7C THt PRODUCT OF 
OlMrNbXUN PX(49C0I 

C THt* FULLO«iM«» OtMtNalON MUST »■» OPEAlER THAN CM tCUAL TC (M^l)#M/2 
OlMt.KSlUN R( 2:00) 

C Thc EULLOWINO C!*»CNS!UN MLS7 gP GPtATER THAN OP EGHAL TC lO 

OtMEKTTON ANSCIO) 

OaMENS t(;N VV ( 7 3) , v*"! 70) 

COMMON/ULNC/CiAi «U • GLPUM • CLPQS • X NOE X 
COMMQN/BLNh/ VV » Vf tANS*TSAV”#M,K 
* FLh»mAT (AA,A2,I5«2I2) 

a MURM AT < l‘^l * 'MULTI (H.E ct<iC-£ESlCN**5X.AA,A2//ex,'£5LECTlLN*,I2.//> 

3 FOFAAT (lM3**VAPIAULE*i‘SX,«M£AN',»>X*»STANDAPD**AX|*CrPPELATlCN** 

I •*K, •FCCP._SbICN**A»**bTO. PPOOO* . jX, 'CCMPOTEC' */*• NL . • * I 6 X , 

■ “ 2 *DEVIATION» .7X,« » vS. v • , 7 X • • CCr:rr I C I ^Nt • * 3x , • OF R C*" . CCFF , • , 3X • • T 

JVALJC ) 

A FOR'ilAT (IH •X<*«bFl4«*) 

5 FORMAT ( XM • 'CcPENOt NT • ) 

3 FJvJMaT ( 1 M^,/,IM , • I NTf RCtPT* , 1 OX, Fib*f • // • Ih .•MULTXPCc CORNCLATI 

ll»N« tFi^. 6 ,// ,1M ,*STO. CPMOM OF T I M ATE • , F I J • 6 # // ) 

7 rOPMA- ( IMC • ^ lx , • ANAL7SX S OF variance FHP Thfc ot C-RF « 5 1 C N • , // , 5 X , 

X 'SJUlCt OF VABlATIf.K*,7x,'DCr,crr^,*,7X,* SUM CF * , 10 < • * M t AN • , I 2 X , 

: »r VALt»F* • 3 3X, • cr f REECOm • , a X , • 33LARE S» , 9X , • SQOAC £ i • ) 
n FC^’•A■' (IM • • A^TPl nU^AMLC TQ RCOf- E S S X ON' , 4 X , X 5 • IF 1 « ♦ f , / , I H , 

I 'Dt_ViAT2JO FRCM P EGRESS t ON ' , A < , 16 , 2Fl t • 5 ) 

^ format (IH , Sx, ' TO TAL • , 1 9X , X o,r I fc. 5) 

10 format (Jcl 2 ) 

11 format < 1 H ,l*;X*'TAnLE OF R€ S T CJ Al. 9 • , / / , 1 H ,'CASE NC.'#''X,'V VALUf 
1 • , UX , ' r c. »TI MA Tt: • , »,x • • RtSl OUAl • ) 

i<f Format c im , ,f i;>, « , 2 Fi a •:> ) 

13' rOHMA"^ (lMl,'KjMb£& OF 3£L?CTlCNr NOT riPECXFItO. JCU Tf Lm T nati n ' ) 
l«* FuNMAT ( 1 ‘^Cf'TOt MATMix IS SINOULAC, ThXS SELECTION '1 S^IPOcO*') 

15 cukMAT (12X,7F13«3) 

C RcAJ THt FPf ULLM PAPAMCTcR CAL'> 

c 

C - ALPHANUMtAlC NAMC jf; •AA‘*‘ELL!T£ 

C N • N'Jmu£P CF OBS£RVArirjfi” 

C M - NUMBER OF VAR IaPLES 

C N 3 - NiJWOtP cr SElTCTIOnA 

c 

UO FoAC ( 3* 1 .LOD- 130)PR .PRl .N.M.NS 
10*0 
X* 3. 3 

MV* M - I 

Ju 30 J-=l*MV 

J 3 VV ( J4* ) 3(. « 3 

CTCLL CrFR” (N,M, I n , X,XrAR , “TR , Rx ,p ,0* V,T ) 
r TFS^T NUMt3tP OF SEL 2 CTXCN; 

IF IN'.) lC4.lC-l.l0i 

136 Ah. lie (o.ij) 

rc *c3 

I 0-* 'JO 200 1 * I .NS 

AS ITE ( >. ^ ) PF .PRi , X 
C RiAO aCBSFT FLLtCTX.lN CA^O 
C ' ■ 

C N - r.UMOCf? CF INOEPLNOENT VAPtAPLES INCLUDEC 

C laAVE - A ViCTjw CCNTAlNlN-'. Tf«c INDEPENDENT VARlABLEr INlLUOcO 
C 

RFAO (j,l7) K • ( I lA VL ( J) • J^l *N) 

C NOeP - l/rPSNOFNT VAFJARlE 
ND£P*«. 9 

CALL t CttP< M.f. ,Nut n, <, I AV=.,RX ,P V) 

CAIL MTMV (R X,K ,Or T. V. T ) 

C TtSr alNGI.'L AR I Tv pF MATRIX iNVERTCn 

Ir ( DuT ) 1 Iw • 1 10 • 1 1 > 

U 3 AM 1 TL ( ' . I < » 

GO TO ^^0 

kid Call mijuTR ( <•<• xd AP ,STu, O.PX *FV , !«;aV£, V ,StJ .T • ANS ) 

C PRINI' M_ANi, >TAN"AAD UtVlATJCNS* ( N T " RC OPP fc L A^ I C NS SCTwEdN X AND V. 

C RFoNt^'lCN LiJiFF 1^. ’ENT , 3TAN0APO C■:VIATTr;^r q|. RcriRtSbtrN CQc FF 1 CX ENT « , AND 
C Cr.MPUTEO T-VAl‘JES 

VV=K ♦ i 
• r 1*^0 U . J) 
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I SN 0047 00 lit J» I .K 

ISN 0048 L«1SAVE(JI 

... ^ S4Vt ntO«4t:S8|ON COEFFICieNTS 4NO THEIR ERRORS FOR SUFRCUTINf LCVNUF 

ISN CC4« VVILI«V(JI 

ISM 0C5C VE(L)«S(MJ) 

iSN OOM - 115 WTm6.41 L.WVreiDtSTDTLl.PVfjl.VfJl.SRf Jl.T< Jl 


X SN 
XSN 
ISM 

0C52 
OCS J 
0064 

c 

4»iTE C6f4| 
L*XSAV£(NM| 
•NITE (6*4) 

P«iNT iNTfRCCPT 

ISN 

0055 


T t **4 T t 

•rite <6*6) 



c 

PRINT 4N4L YSI S ; 

1 



•W XT C TC • T^l 

ISN 

C057 


U«AnS( 6) 

ISN 

O^Sfi 


4kIT£ 

t SN 

0C5<; 


U«N - X 

ISN 

3060 


SUM«ANS(4) 

t SN 

0061 


•RITE lOfOl 

T SN 

0062 

c 

CALL bUURUUTlNf. 
CALL LOVNUM 

-Tsrr 

-peg"? — 

~ ~2in — covriFnjT" 

I £N 

CO^ 4 


GC TO 100 

ISN 

0065 


300 CONTINUE 

ISN 

CC66 


h£T<JRN 

!£N 

0067 


tNO 


► 4NS(7) 

L I SOM 


62 


I 


ItN COOi 


'>.J[iPC>liT IN« 7»-FW3V (JPLOT*»^^l 


1£N Cr3c 
1 3N OoO^ 


! SN 300= 


I jN 

C3 7F. 

I SN 

3 337 

I ?N 

003“ 

I iN 

000-^ 

I iN 

c <: 1 0 

♦ iF' 

0 31 1 

ION 

371 ? 

I £N 

3013 

: SN 


t On 

001 s 

1 fN 

CGI c 

: =N 

COl 7 

fSN 

031 S 

f ?N 

C j 1 ■? 

1 SN 

0C2 C 

I ,»N 

oo;ii 

? £N 

0322 

1 3N 

G02 3 

' TSrr 
I SN 

0'*24 

002? 

liN 

0 32 6 

? ^N 

C 027 

I SN 

0 C 2 6 

I af. 

0029 

1 FN 

0C?0 

- -TTN-CTT * 

IFN 

00 JS 

I £N 

0333 

ICN 

C03^ 

I £iN 

003* 

TtN 

xmr 

I SN 

0 33® 

I IN 

0 C ^ 0 

ISN 

OOA c 

I SN 

00<«1 

ISN 


1 SN 

00 43 

I SN 

0 C 44 

1 TN J ’J«*0 

I TN 

0C46 

ISfJ 

0 34.- 7 

I ?N 

0 0 V »j 

I FN 

CC49 

IF.N 

OCbO 

I SS 

CCbl 

! £N 

0C52 

-T-tTf -TTtT5T‘ 


,<r ^,n’K^Ttn- cj"PUTLS thf •mplit’jct'- r- tmf ccnktiitlfnt 

IP th' T *P r in Th- 'OlID patTh, GIV^N THt fCLiu) Fac.m i. ivE N*J^PF<5# 
ort ii lo-H^ zl”. H.Jtl !T ALFO plots -^..OLTINC TJP AS A 


TMI 

T1 

LA»i AlvOL t iS _ - - 
jAmOc- PL.-T ir De*"!'-*!.? (JPL0T = 1 


If vr« 


ijcTc. LtiwiTOLAP 


5 Kl ANP f'2 A^e L'JMPcO IN WiTH ‘^f'* LJ‘MQ 


NOTAT 'UN 
P u HD^TL'^tLPtLQ) 
( uM*Lp.i-i ) 
AMjJML^iLPiUy) 
A* t L'» tL o » 


PC^nD OF LUNAP CTKjTITlIENT LW,lP*LG 

- period JF -iLLAR C3N3TITJENT L»'*LP.LC . o /• e r-/*r .. r.e 

- AMPLITUOL or LUNAR CONSTITLtFT lv.lP.LO IN ARCSECCNDS 

AMPL ITUOL OF 'OLAR CCNLTITUrMT LW.LO.uO IN ARCSECrNL-S 


I.ARaUF(2f0.-l) tAR 


jMPi.lCP' RE A L* “ < A-“f*t ^ 3 
Jl ir N i.PSL»:M ( 2 . 1 * *» ) • KPr^LN^i ( i 

1 APGCTTMC 3 • 3 • 3 3 t Arr,cT;( 2* 3 • T) i s i i » 

-IM^Ni ICN LntE».««OCS(2. S.?) 


, 4»-< . 
, t»-K • 


I Mu s' 4 T? «p ( JO J ) , xie<P( ?03> • T (7 00 ?• XI NCL ( 7 00 » •N’^K AC*-. ( 7 00 ) 

I-JN C( 2, 1) ,AI (2, ,U<2.3).RLPCV(3. M .GLPO*- (i, 3 > 

DIMcK. ICN TI or ( 2.3* 3 » 

kiiMt.NiK*N :JAT(1 h» 

DiMuSElDN DuRPM( 2, J, *» I , ( i 1 1 • 3 ) 

r I5MM "'^’/t■^LK A/ F F >L ^ P«- L N H • AP GU*^ • AP G'J C • AP GO fM t AR GL 

CTM-»5*.Vnt;KB^lj:^W.^FNTii.CFNTS.SU‘'.XMC0N.FSUN,TFN.VFS.VFt.- 

L"v:5;:;;Li:Ju;iLr.;u'n;r:i;:iF^u?J«^^^^ .<ic.A.L.Ki.xNcn€ 

l T 7 ERfl.Kl MLOS . X! SUN. TMS . XNOM ,y fjpt t F. XMr. AN'* • X M£ ANc. X C' T‘*, XMOn , 
c GMGOLT,M£xP 

uCMMCN/tiLKG/C # Ai tn*(»LRQW,bLPOG.lS3£X 
CL' MMON/HLN I/S*- 

DATA TM II 2 CTMI ’2 .TM 2 12 .7M2 ->?/AMO 1 t AHK I ,4f*M2 

L.ATA 73112. F3l2i .T52l2tTF2?2/AHni , amf 1 3 .#f*S2 

C PlNij* 7 mJ'i Ai^^“ "7 'A^PL?‘*’'IOt tJF ALU ’Hf CONS'*’ I TUF NT S (IG^CKI^£ ANY CCNS"^ I'^NtST 

c -i?« i PSP ?-o:nNOuR than ?300 lays, the ?000 is an APR. Th ary FIGU..E.) 

AMAX^C «0 
0 U 'j 3 L M = i , ? 

YL** = LM - 

OCT- CO 

JU ^3 LO=l . s 
YKM( L^fUP* *LO ) =xs 

XMT OS t LM.LP.LC )=YKV(Lf4.UP.LC)»CFSTM«ni(UM)'C(LM.«:)*Al(LY,LP)*‘ 

1 N»Lf* C'-(LPiL‘i)*XLM/(AOGDTM(LA*LP.LO)*SS> 

X M COS ( lM, LF, L 0 ) =DAUS< XMCCS ( L V . LP .LQ ) ) , 

P£RDM(LM.lP.UO|- 2.0‘PI /( ARGr‘'’V(L‘'.LP.LO» » 

-T»p-mw tt-v. 1 PrT_Tti sDARS(PEPOM(L**,LP .UO) ) 

Tf^U^LP.^^)=YXsVL^ .CFNTG.RULM)-C(LM.2 )«c(LM.lP|. 

1 GLPCL(LP.L0)*TLM/( aqgO^S(LM,LP*L ) J«S5) 

X^UN(LM,LP*LC)=CAUS( XCUN(LM,LP.L0I ) 

c e> x«ooN(L«.LP.LQ)=«MorNUY.LP.La. « 

1 X£'IN<LM,LP.LG) ^ . 

P£R05( L^-. LP.LQ ) = 2 .)»Pl/(AfiGr*T->(LK.LP.LQ)> 

t»FROTrLV.l-Pt LO) -0 A0S( PEPOG ( t. v.LP tLOl 1 

Pt.M=PLRf. M(LM«L°.LG) 

IF (P^R-JJOO.O) 9**ff.4,62 
c4 COM INUE ^ ^ ^ ^ , 

If < XMCCNC LM ,LP .LQ )-XMAX ) 62,6«L.t3 

0 i xmax=xmooMlx.lp.lo) 

62 CuNTlNUL 

PfcP=PrPOG(LM .Ln^.LO) 

tF fPTP - 2 C 03 *01 *H.C5.t6 

6 5 CUNTINIJL 

IF ( Xf'JNlLM. LP.LO) - XMAX) 60 . Co.'S? 

07 XMAX=XbUN(LM .L^.LO) 

©6 CONTIMI'L 

60 CDNTlNUf 

ITS < C . J ) 

1 format { lHi.///> - 

'.^oiTr'te • I ) 

c pSINT^)JT''L^^Ip°ot C'JNSTITijlNT N AME . DI 5TURO I NG QCDY. AYPLITUGE IN ARCSEC. 
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C W^aiul) IN CAri , OiU CfNT or CAaoEj-T CON^TITUSNT, and !fAV::(J» VALUES 


I SN 


1 

rr;f;wAT (• i. n i> 0 t-ql niP.ufoY amplitujl 

PERIOD 




1 PcR Cf.NT ISAVclJ) (FQP 'IBF ?N SUrlPCUTlNt PLGPESIM 


: 3N 

0 


1 f e ( C, i ) 


! £N 

0 0£>e 

2 

Format ( 32*# #10X#» (CAV*) •,/) 


ISK 

rcvr 


on 39 LM^l ,2 


I IN 

CCS® 


DU 39 LP=l»' 


: £N 

0C^9 


UU 09 LO= 1 #3 


T fN 

C C 0 

0 4 

TIDlILY.LPiLCI =IjLANA 


I 5N 

OOsl 


TiDu u • i . : 1= rvi i#s 


i ?N 

00*? 


TiDe ( 1 t 2# 2 )eTMl >2 


’ jK 

C Cf- 3 


TID-I 2. 1 #2 ) = TM21i. 


T CN 

OCt 4 


TIC-t(rt^#2 »*TMj2i 


ITK- 

rns* 


LLL^2 


1 SN 

C Coc 


1 = 0 


I SN 

0 0c*? 


DO 70 L‘1=1 ,2 


T ‘‘N 

0''t = 


OC ’ J LP=1 . 1 


r 

oc* •; 


LPP=LP - 1 


1 £>N 

C07C 


OO *>0 LO = l #5 


1 SN 

CC71 


LOU=LU - w 


! cN 

00-2 


AMPS (XMC;ON(L»'fLP#LO) /UMAX! ^0.0 



0C73 


I=I *■ 1 


! f.N 

0C7<. 


I 3=1 


T e N 

0 :7i 


1=1*1 


I a>N 

0 07i- 


IF (LPI* .LQ* J •AND# L.O'J .SO* 01 “»C *^0 71 


* 5 N 

CC7t 


aPITL (jtj) LLL .LM iLPP.LCC .T IPF. (LM ,LP ,LO 1 , XMCOM LM .LPiLO) . 





1 PERa^'lLMi LP *LO» ■ AMP# 13,1 


IZr 

C C '**'* 

l> 

rORM AT j, ?*,A4,FX#* MOCN* #'iX*JlT.D#JX#FlO#3#3X#Fl0.1#4X# 

12# 




1 • AM • # i 2 » 


ts:m’ 

OCPC 


Gij ro ’c 


I SN 

C Of J 

71 

YPITt (c.'.) LLL iL M .LPP.LQQ .Y 1 DE ( LY .LP .LO » , X VCOMLM.LP tLC) , 





1 PcPOML.y# LC* #L0J #AMP,1 3# I 


! £N 

V 

o 

F (4 I c, 2X , Aa, 3X . • MCnN*EUN '#3X,Dl l#A,-3x#F10*l#’X#FlC.3 

.UX . 13 • 




\ • ANC • . I 2 ) 


! fN 

0C62 

70 

CONT !^■Ur. 


1 CN 

COBh 


«P ITF ( 0# 4 ) 


I £N 

OC BS 

4 

FORMAT ( 1 OX# /✓ ) 




C JO^HE Tire? 




C PhlNf uljr L.M.P.G. CDNSTiruEMT >1 AMc , CI STL'PH J NG DCOY, AVFLITUCE 

IN ABC5SC* 



C f'tPlJD IN CAYS. PEP CENT DF LAPGEST CONoTITUcNT. AND lEAVc(J) 

VALUES 

I SN 

CCBc 


on ie LMsi 


I SN 



DO -36 LPsl . 3 


f SN 

o>-^« 


DO -3® LO = l #3 


I iN 

0CS9 

o« 

T IOC (L M#LP #L r.) =«CAN« 


I iN 

OOc-0 


TU'e < 1 # i # 2 )= T£1 1 J 


1 hN 

7^/1 


T I DF- (i#.j#.?)=Tbl22 


I SN 

0092 


Tine ( :# i #2 ) = 1 S2i2 


I^N 

C 0^3 


TIL'c le#2#2 > = TS2i:2 


I SN 

c: 4-* 


on tio c^ = i # 7 


I 5N 

oc**: 


oa 4c LP=i . j 


! SN 

V V 


LPP=CP - 1 


: SN 

c:o7 


or eo Lo =: #2 


T SN 

oc=«e 


LOoscQ - 2 


I SN • 

C 0 IV 


ATnT=C3T?yN(LM #LP#LO)/XMAiO*lCO.O 


I SN 

OiOO 


IF (LPP •£'3# 1 .ANO# LOQ .EQ# C) '30 TO A1 


I iN 

0102 


I = I ♦ 1 


I £N 

CIO J 


I 3=1 


! 5N 

0104 


I=! *1 


I SN 

OIOS 


• Pile (t.7) LLLiLM.LPO.LOO.T I Oe<LW iLP.LO) .XSUMLY.LP .L')» t 





1 Pc-POC (LM# LP #LOI ,AMP , I 3, X 


ISN 

010c 

7 

rjKJAT (Al2,3x,A<,5x,«EUN*,f, X,Oij.6,3X.F10.3.3X,FlC.;.AX.I2, 




1 • AMP • • 12) 


: SN 

c i : 7 


GC TO 30 


!SN 

oioa 

01 

XPITE (D.3) LLLiLM.LPP.LOQ.TlDECLM.LPtLQ) 


I 3N 

0 109 

*4 

FDPMAT < AI i, 3X. A4, tx , • SUN* .1CX.*.,.>.1AX.*...',11X.'...',‘S 

X# • • • • • ) 

I £N 

Ci 1 C 

00 

CDNT INUr 







i SN 

C 1 I £ 

?0 

FOPMAl ( /////.lOX, •TH cORETICAL amplitudes assuming bCLiD EARTH LOV 




IE N'JMCeF K2 = • #F7#3 1 


1 SN 

U 1 1 J 


•WW I Tt TT., 71 ) COATfJ)# J = 1#IA) 


1 SN 

0 114 

21 

FOPMAT </////#lOX#13A7#A2) 




C PWlNT CUT A PLUT CF THE SCJLlD -ARTM TIP t?F JPLCIxll 


ISN 

oil s 


IF IJPLCr #ec. 1) GC TO 63 


1 SN 

0117 


GL TU 12 


I iN 

OllP 

03 

CoNTiNUfc 


I SN 

0 1 1 V 


JG 10 L = 1 # VT 


ISN 

0 12C 


TT=T(C ) 


T5FT 

oizr 


CALf Al3f5(TTl 


I SN 

0 1 ?? 


SUM=0.0 


ISN 

.0121 


OU 90 1_M=1 #2 
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lEN 

0 1?A 


YLM=LM 

I SN 

0125 


00 90 LPsl • 3 

XSN 

012F 


UU 90 LO«l *3 

I5N 

Oli 7 


XPC‘ON(LM.Lr. LQ) =YK M ( L M , LP . LO ) «CFnTm* 01 (L^>»C(LPf2)*Al ( L P . L P ) »GL POP 




1 (LP« LU) ‘YLP*^ AP SNM(LM,LP.LQ ) / ( ARG P ( L M , L P . L 0 )«5S ) 

ISN 

CT^3 

' ~ 

3n:UVTi:iif.LP .L CI^YKS (lP.LP .LO) fcCFMTS«rTl(LM> kCUP.2)«E(Lik.LP) •glposcl 




IP.LQ)' YLM* ARGLS(LM,LP*LU)/( APGCTS(LM,LP#L0)*SS) 

ISN 

0129 


dJHsaUM * XPOCK(LMaLP.LQ) f XSLN (LP«LP .LO) 

I 5N 

01 (C 


90 CONTlNUf 

I SN 

0131 


X1NCL(L)=SLM 

I SN 

01 1c 


10 CtiNTlNUE 

2 SN 

0 1 33 


call PLCiTER (MT.PcXP) 

2 =N 

013A 


12 CONTINUE 

I SN 

OliE ■ 


RFTUnN ■ ■ 

I SN 

0136 


cNO 


1 SN 


! ?N 0003 
I tN occ*» 

I 

ihi c c 3 7 

iSN o^ce 

2SN OCO^ 

I £N :oio 

OOil 


bOUPOUTiNr. T t^F t (TSTAHT, TCKO *rT * •*) 

C 

c 

C ThI5“*5TjHK^UT IK? FILL? IN The VALUES OF A«PAY T<L). L»V KNOWING ThE STAPTING 

C TiHc (T«.TAPT), TMf ENDING TIME (TENT), AND THE TIME INTERVAL t*ET»EFN 
C PJINTS (CT). 

C 

C 

implicit PEAL*« <A-M*n-»’) 

01 PEN* r f)N TF tPIi'OO ) • »I EXPC »C0) .T (■'C'? ) • XI NCL ( 7 00 ) .NT S AC" (700) 

CJMMGN/nL«.C/ FE XP tXlC XO Xl NCL .N^PACK 
M=(TcND-r>TACT )/0T ■*• I.IDO 

OL 12 L = l • M 
XMsL-i 

12 T(L)=TSTAPT «. XM*[)T 
Pt TURN 
CND 
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WCTN 
^ UN 


aTcLLITu -icTi - I t l5'o-l5hT cata 


satellite L'ATA 


A 

{ 1 p C-1 ) 

£ 

INCL INAT JHM 

roESPFE*:) 

NOCE PATF 
(OEG/CAY ) 

T7F«a 
( CAYS > 

NODE 

CCEGREE3) 



0. 07<ii. 

55. 3fl05 
MCC^ ANe"> 

-2.?4< 5 
3UN CATA 

36.5000 

-I0<5.0fl20 


OMEGA 

OLLTA 

! NCL 1 NA*^T CM 

TM< 

MEAN ANPMALY 
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Table 1 


Orbital Data for the BE-C, GEOS-I, and 
GEOS-Il Satellites 


Satellite 

Semimajor Axis 
(10“ cm) 

Inclination 

(degrees) 

Eccentricity 

Nodal Rate 

0 

(degrees/day) 

BLi-C' 

7.5o:: 

41.1667 

0.025 

4.2535 

GEOS-1 

8.0729 

59.3805 

0.073 

-2.2465 

GKOS-II 

7.7052 

105.7896 

0.033 

1.3997 
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Table 2 


Tidal Parameters Derived from the Multiple Linear Regression Analyses of the Tidal Perturbations 
in the Orbital Inclinations of the BE-C, GEOS-I, and GEOS-II Satellites 
(The errors indicated are the formal statistical errors. ) 


Satellite Tide Cnipq 


C.HOS-I 


Cnipq 


Period 

(days) 

Amplitude 
(arc sec) 

^2mp0 

*2mpO 

^2nn'0 

^ 2nipfl 
(deg) 

(cm) 

2100 


11.8 

0.09 

0.276 1 0.077 

-0.018 10.088 

0.28 • 0.08 

- 3.6 1 18.1 

3.0 1 8.1 

2110 


84.8 

0.80 

0.254 1 0.010 

0.025 1 0.007 

0.26 ±0.01 

5.7 1 1.6 

6.0 1 1 .2 

2200 


10.3 

0.11 

0.204 1 0.053 

0.030 1 0.054 

0.21 ±0.05 

8.4 1 15.1 

10.3 1 5.4 

2210 


42.4 

0.12 

0.322 1 0.065 

0.128 1 0.072 

0.35 X 0.06 

21.7 1 1 1.6 

2.1 1 1.2 

2100 


57.8 

0.13 

0.205 1 0.036 

-0.016 ±0.041 

0.2 1 1 0.04 

- 4.5 1 1 1 .4 

3.8 1 1.5 

2200 


34.4 

0.20 

0.253 1 0.030 

0.025 1 0.030 

0.25 1 0.03 

5.6 1 6.8 

2.5 1 1.4 



0.241 10.005 


0.251 ±0.027 


0.284 ± 0.034 


0.181 ±0.044 


0.2^7 ±0.01 


0.011 ±0.005 0.24 ±0.01 


0.017 ±0.027 0.25 ±0.03 


0.000 ± 0.035 0.28 ± 0.03 


0.010 ±0.045 0.18 ±0.04 


0.034 ±0.01 2 1 0..30±0.0I 


.7 ± 1.3 7.9 10.7 


3.9 ± 6.1 5.3 12.8 


0.0 ± 7.1 0.3 10.6 


3.3 ± 14.2 4.8 ± 1.8 


1 .7 ± 0.6 


0.257 10.021 -0.046 10.013 0.26 ± 0.02 -10.2 ± 2.7 84 1 2.2 

0.242 1 0.025 0.050 ± 0.025 0.25 ± 0.02 H.7± 5.7 7.8 ± 2.5 

0.306 1 0.023 0.012 1 0.029 0.3110.02 2.2 ± 5.3 0.2 ± 0.5 

0.050 1 0.022 -0.107 1 0.037 0.12 ±0.03 -65.0 ±12.3 10.9 11.0 


0.340 10.007 -0.015 10.004 0.34 10.01 - 2.5 ± 0.6 2.110.3 


(deg) 
144 1 163 
29 1 8 
287 1 31 
101 28 
1 70 1 24 
2‘)8 1 33 


1115 


289 1 29 


185 1 22 


355 1 20 



313 1 16 


8 


8 1 105 


157 1 7 


11116 


























































































Table 3 


Secular Changes in tne Moon's Orbit and the Earth's Orientation 
and Spin, and the Dissipation of Energy for Selected Tides 
(The entries are a weighted average over all three satellites 
The errors indicated are the formal statistical errors ) 



2100 

2110 

2200 

2210 

m 

2110 

2200 

2210 

"m 

(arc scc/( I00yr)“) 

3 + 13 

0 

-29 ± 15 

0 

- 

- 

- 

- 

"m 

(cin/yr) 

-0.4 1 .9 

0 

4.3 ± 2.2 

0 

- 

- 

- 

- 

^’m 

(/I0‘' yr) 

0.0 + 0.001 

0 

-0.002 ±0.001 

0 


- 

- 

- 

j 

(deg/ 10^ yr) 

-0.3 + 1.5 

-0.2 ±0.1 

-0.7 ± 0.3 

-0.0 ±0.1 

- 


- 

- 

i 

s 

(deg/K)*’ yr) 

-1.9 ± 9.2 

-0.7 ± 0.3 

3.2 ± 1.6 

-0.1 ±0.2 

-1.0 ±0.5 

-0.2 ±0.1 

-0.2 ±0.1 

-0.0 ±0.1 

0 

(10'^^ rad/sec* ) 

0.3 + 1 .4 

-0.1 ±0.1 

-6.3 ± 3.2 

-0.0 ±0.1 

0.1 ±0.1 

-0.0 ±0.1 

+0.4 ±0.2 

-0.0 ±0.1 

i; 

( 10' erg/scc) 

0.1 ±0.7 

-0.1 ±0.1 

-3.6 ± 1 .8 

-0.0 ±0.1 

0.1 ±0.1 

-0.0 ± 0.1 

+0.2 ±0.1 

-0.0 ±0.1 








Table 4 


Secular Changes in the Moon’s Orbit and the Earth's Orientation 
and Spin, and the Dissipation of Energy for Selected Tides 
(The entries are a weighted average over the BE-C and GEOS-I satellites 
The errors indicated are the formal statistical errors) 





(Moon) 



ixiw. 

(Sun) 


Knipq 









X 

2100 

2110 

2200 

2210 

2100 

2110 

2200 

2210 

"m 

(arc sec/(100yr)") 

' 

3± 13 

0 

-17 ± 20 

0 

- 

- 

- 

- 

a. 

(cni/yr) 

-0.4 + 1 .0 

0 

2.4 ± 3.0 

0 

- 

- 

- 

- 

s, 

(/lO** yr) 

0.0 + 0.001 

0 

-0.001 ±0.001 

0 

- 

- 

- 

- 

j 

(dcg/IO" yr) 

-0.3 + 1 ..S 

-0.3 ±0.1 

-0.4 ± 0.5 

-0.0 ± 0. 1 

- 

- 

- 

- 

i 

s 

i(lcg/I0‘' yr) 

-l.‘) ±‘).2 

-1.3 ±0.3 

1.8 ± 2.3 

-0.2 + 0.2 

-0.1 ±0.7 

-0.3 ±0.1 

0.7 ±0.2 

-0.0 ±0.1 

0 

(l(r^’ rad/scc^) 

0.3 + 1 .4 

-0.2 ±0.1 

-3.6 ± 4.4 

-0.0 ±0.1 

0.0 ±0.1 

-0.0 ±0.1 

-1.3 ±0.4 

-0.0 ±0.1 

i; 

(!{)''’ erp/scc) 


0.1 +0.7 

-0.1 +0.1 

-2.0 ± 2.5 

-0.0 ±0.1 

0.0 ±0.1 

-0.0 ±0.1 

-0.7 ± 0.3 

-0.0 ±0.1 



Table 5 


P for t = 2, q*= 0 for the Quantities Shown in the Left-Hand Column 

Cmi'ii 


’’Cmnu 



2100 

2110 

2120 

2200 

2210 

2220 

2100 

21 10 

2120 

2200 

2210 

2220 














"m 

(arc sec/( 100 yr)^ ) 

-144.6 

0 

+0..^ 

-847.1 

0 

+0.0 







(cin/yr) 

+ 21.4 

0 

-0.0 

+ 125.1 

0 

-0.0 

- 

- 

- 

- 

— 

— 

(/10‘* yr) 

- 0.0076 

0 

+0.0 

- 0.0447 

0 

+0.0 

- 

- 

- 

— 

— 

— 

j 

(dcs/lO" yr) 

+ 16.7 

-18.2 

-0.1 

- 1‘).4 

-1.9 

-0.0 

- 

— 

— 

— 

— 


i 

(cleg/ 1 O'* yr) 

+ 10.‘I..S 

-8 1 .0 

-0.6 

+ 94.0 

-8.4 

-0.0 

+22.5 

-17.5 

-0.1 

+ 19.9 

-1.6 

-0.0 

0 

(10‘ ‘‘ rad/scc‘ ) 

- 1.S.6 

-14.2 

-0.0 

-182.8 

-1.5 

0.0 

- 3.3 

- 3.1 

-0.0 

-38.7 

-0.3 

-0.0 

\ 

( 10‘‘' erg/sec) 

- S.-S 

- 8..^ 

4).0 

-103.6 


-0.0 

- 2.0 

- 1.8 

-0.0 

-22.8 

-0.2 

-0.0 






RESIDUALS (ARC SECONDS) 
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TIME lJUUAN DATE! 


Inclination Data and legression Curve for GEOS-I 







TIME WUtlAN DATE) 


Figiu’e 4. Eesiduals for Regression Curve Shown in Figure 3. 
Standard Error of Estimate: 0. 132 
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Figure 5. Inclination Data £ 
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Figure 7. ^herieal Triangle Used for Compiiting the 
Moon's Position With. Respect to the Eax-th's Equator 


98 


